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Abstract 

Accretion processes in quasars and active galactic nuclei are still poorly under- 
stood, especially as far as the connection between observed spectral properties and 
physical parameters is concerned. Quasars show an additional degree of complexity 
compared to stars that is related to anisotropic emission/obscuration influencing the 
observed properties in most spectral ranges. This complicating factor has hampered 
efforts to define the equivalent of an Hertzsprung-Russel diagram for quasars. Even 
if it has recently become possible to estimate black hole mass and Eddington ratio 
for sources using optical and UV broad emission lines, the results are still plagued by 
large uncertainties. Nevertheless, robust trends are emerging from multivariate analy- 
sis of large spectral datasets of quasars. A firm observational basis is being laid out by 
accurate measurements of broad emission line properties especially when the source 
rest-frame is known. We consider the most widely discussed correlations (i.e. the 
so-called "eigenvector 1 parameter space" and the "Baldwin effect") and analyze how 
they can be explained in terms of accretion properties, broad line region structure, and 
source evolution. We critically review recent estimates of black hole mass, accretion 
rate, spin and possible orientation indicators, stressing that any improvement in these 
parameters will provide a much better understanding of the physics and dynamics of 
the region producing the optical and UV broad emission lines. More accurate mea- 
surements of Eddington ratio and black hole mass may have a significant impact on 
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our ideas about evolution of quasar properties with redshift and luminosity as well as 
on broader cosmological issues. 

1 Introduction 

Tremendous data gathering advancements make now possible to see a distant solution for 
most of the deepest conundrums concerning quasar research. There has been an exponen- 
tial growth in the number of known quasars lasting since the 1970s. Presently, the 11*^ 
edition of the Veron-Cetty & Veron Catalogue of AGNs and quasars 12531 lists «50000 
quasars. The Sloan Digital Sky Survey (SDSS) plans to compile a sample of ~ 100000 
quasars. The SDSS Third Data Release provides Charge Coupled Device (CCD) spectra 
for 51000 quasars found over 4200 deg^ 12101 . The improvement in observational capa- 
bilities has been possible by the introduction and spread of CCDs and other linear, high 
detective quantum efficiency devices as detectors for astronomical observations since the 
late 1980s. At a second stance comes the increase in access to large light-gathering power 
telescopes of aperture ^ 4m, crucial for spectroscopic observations. The ability to carry 
out multi-frequency observations with increasing spectral resolution and sensitivity (most 
notably provided by Hubble Space Telescope (HST) and by the Far Ultraviolet Spectro- 
scopic Explorer (FUSE) for the optical/UV) through the 1990s and early 2000s has been a 
third factor of relevance. The data gathering improvements - which are the foundation of 
every astronomical advancement - have led to the discovery of systematic trends in quasar 
properties. These suggest that quasars are not well described by an average spectrum. We 
can amplify this statement to say that the spectra of quasars seen at a fixed viewing an- 
gle are also not the same - the basic tenet of Unification Schemes 1 1 1. Yet, data have not 
been fully digested in physical terms. Observational constraints have been only very re- 
cently organized in meaningful ways. We still lack the ability to derive important physical 
information from observational parameters on an object-by-object basis. The first aim of 
this paper is to set the point of our present knowledge on the diversity of quasar spectral 
properties, especially in the optical and UV spectral ranges (§ We then discuss much 
needed improvement to ensure accurate measurements of the main physical parameters (§ 
EJl as well as structural constraints on the line emitting regions (§ ISj and the physical basis 
of quasar diversity (§|5j §0 §[8ll. 

1.1 Basic Accretion Parameters 

We regard quasars and active galactic nuclei as systems accreting matter onto a massive 
( ^ 10^ Mq) black hole (see ^1.31 for possible caveats). The main parameters are, as for 
any accreting system, black hole mass (Mbh» ®. Eddington ratio (©, and mass accretion 
rate M in Mq yr^^. The Eddington ratio Lboi/^Edd is defined as the ratio between the 
AGN bolometric luminosity Lboi and the Eddington luminosity L-^dA, i-e-> the limiting 
luminosity beyond which radiation pressure overcomes gravitational attraction if accretion 
is spherical II177I . The Eddington luminosity is directly proportional to Mbh> and can be 
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written as 

LEdd ~ 1-3 X 1O38(Mbh/M0) ergs s-\ 

Since the power emitted by an active nucleus through conversion of mass into energy 
can be expressed as -Lboi= "qMc^, it is possible to define an Eddington accretion rate 
-^Edd = -^^Edd/(??c^)> and hence a dimensionless accretion rate fn = M/MEdd- We 
stress that LboiZ-^^EddOc LboiA^BH> and that Lboi/AfeH is a quantity that can be derived 
from observations. On the other hand, the value of the efficiency 7] (~ 0.1) depends on 
the accretion mode which may not be the same for all sources: if matter is confined in an 
accretion disk as it is customarily assumed [211], rj depends on the geometry and radiative 
properties of the disk, which in turn may depend on M). Therefore we will not always 
confuse the dimensionless accretion rate m and the Eddington ratio as it is frequently done 
in literature: for a fixed efficiency, there could be well super-Eddington accretion even if 
the source is radiating at, or below, the Eddington limit |41 1. Black hole spin and an 
orientation angle (defined as the angle between the line of sight and the axis of the accretion 
disk around the black hole 0, ^5.2t probably matter in the context of UV/optical properties 
although they have turned out to be very elusive to measure. 

1.2 Nomenclature and Samples 

We use the word quasar here in a generic fashion which means all classes of extragalactic 
objects that show broad (full width at half maximum, FWHM > 1000 km s'^) optical 
and UV emission lines. This includes the nuclei of Seyfert and broad-line radio galaxies 
(BLRG) as well as radio-quiet (RQ) and radio-loud (RL) optically unresolved sources. They 
are often referred to under the umbrella of Active Galactic Nuclei (AGNs) which reenforces 
the paradigm that they are driven by the same physics differing only in their redshift-implied 
distances and, hence, luminosity. At the same time, there is no divide in AGN occurrence at 
the canonical boundary (absolute B magnitude Mb ~ 22.5 for Hq 65 km s^^Mpc^^) that 
separates Seyfert nuclei from quasars 12531 . so that nomenclature may be kept luminosity- 
independent. For instance, we keep using the word Narrow Line Seyfert 1 ( ^2.7t for sources 
that are actually luminous quasars. 

Most sources considered in this review were "classically" selected through color crite- 
ria. Major color-based surveys include the Palomar-Green (PG), the Large Bright Quasar 
Survey (LBQS), the ESO-Hamburg (HE) quasar survey and the SDSS. The SDSS photo- 
metric system was designed to allow quasars at ^ z ^ 6 to be identified with multicolor 
selection techniques |208|. Other high-quality data-sets, even if more heterogeneous and 
less complete, are also considered l22l l50l l89l 11401 l2T2l I277II . especially if spectral resolu- 
tion A/AA ~ 1000 and continuum S/N > 20. These samples often include ~ 100 objects 
of "good" spectra; studies on the LBQS involve ~ 10^ sources, while the SDSS accounts 
for another order of magnitude leap in sample size, with ~ 10^ sources. 

No consideration to obscured AGNs will be given since the very possibility of accretion 
parameter estimation from the observed optical/UV quasar spectra turns out to be related 
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to the measurement of quasar broad emission line shifts, profile widths, and equivalent 
widths. Quasars showing broad lines are often referred to as "type-1" AGNs to distinguish 
them from Seyfert 2 nuclei and type 2 quasars (see e.g., Refs. fl1 ll73l ) which will not be 
considered here. 

There aie more cumbersome, borderline AGNs which could be grouped into three 
classes: 

1. sources that show very large broad Balmer intensity decrement Ha/H/? ^ 7 |61|. 
Objects of this kind often show dramatic differences between the broad component 
of Ha and H/3, and are not always found in color-based surveys. RL AGNs hosted in 
early-type galaxies like Pictor A are among them I228II . 

2. An analysis of spectra of quasars in the Half-Jansky Parkes Survey l76l suggests that 
the wide majority (ps 80%) are sources whose continuum is dominated by a pure 
power-law spectrum {fy oc v^~- ^) ascribed to synchrotron radiation. These sources 
show broad lines and are bona-fide type-1 sources although their equivalent width is 
somewhat lower than that expected for color-selected quasars (see also Sources 
with rather low W(H/3bc) and very broad HellA4686 (most notably NGC 1275, in- 
deed a core-dominated BRLG) may be seen as somewhat peculiar with respect to the 
other AGNs, especially if an important goal of the spectral analysis is A/bh measure- 
ments. 

3. A third class of sources may show contamination in the nuclear spectrum of emission 
components associated to strong star formation. At the very least, the bolometric 
luminosity can be significantly affected as in the case of Mkn 231 l26l l53l [2251 . 
Circumnuclear starbursts do occur in Seyfert 1 sources and quasars f93l ll02ll57l . and 
are usually not resolved. The emission line spectra of luminous Seyfert 1 nuclei may 
well be affected by the absorption/emission features observed in strong wind galaxies 
(e.g. NGC 4691; IHTTl ). Broad Absorption Lines (BAL) observed in quasars could be 
associated to nova stars 12191 . 

We mention these three AGN optical/UV spectral typologies since they may be ex- 
tremes of AGNs accretion parameters, especially of .^^bol/^Edd (SI sources in the 
third class are likely to be "young" quasars radiating at high -Lboi/^Edd; the first class 
may include "dying" quasars radiating at very low -Lboi/^Edd with little reservoir of matter 
(gO- 

1.3 Rationale for a Black Hole in Active Galactic Nuclei 

Do we have a definitive proof that the central massive object at the center of quasar is a 
black hole? From the definition of event horizon (i?g= 2GM/c^ for a non-rotating black 
hole, where G is the gravitational constant), the mass of the black hole must satisfy the 
condition M^Y^IRg ~ 6.7 • 10^^ g cm^^. Although this criterion has not been satisfied yet 
by observations of extragalactic sources, circumstantial evidence in favor of a black hole is 
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Figure 1: SDSS median composite spectrum, in the UV (left panel) and in the optical 
spectral range (right panel). Major features are identified; regions with strong emission 
blends of Fell are shaded. For a more exhaustive line list see Ref. 12471. 



now considered overwhelming fT6l. From the time delay in the response of emission lines 
to continuum variations ("reverberation mapping") it is possible to estimate the central mass 
assuming that the line broadening is due to Doppler effect, and to test whether the emitting 
gas motion is virial. In luminous Seyfert 1 galaxies, the central mass can be 10^ Mq 
within the emitting region of the broad lines (the "Broad Line Region", BLR) size, of 
the order of 1 light month. The resulting Mbh/'^blr~ 10^^ g cm^^ is far less than the 
Mbh/?^ requirement for a black hole. However, time responses for lines of different width 
are roughly consistent with a Keplerian trend i.e., with a "point-like" mass concentration 
with respect to the BLR size 1 123 1. It is therefore legitimate to talk at the very least about 
"compact massive objects" at the center of AGNs. This would be a proper, observation- 
bounded definition. We will use the term black hole with the reassurance of a large body 
of cirumstantial evidence as well as with the perspective of decisive tests. In addition, most 
of the reasoning presented in this paper is largely independent on the exact nature of the 
"compact massive object," provided that it is really compact and massive. 



2 Eigenvector 1 

2.1 Quasar Surveys and Quasar Spectral Diversity 

Fig. Qshows a median composite quasar spectrum, between 800 A and 7000 A, with sources 
spanning a redshift range 0.044 < z <4.789. The spectrum was obtained using a homoge- 
neous data set of over 2200 spectra from the SDSS, at a resolution A/ A A « 1800. It reaches 
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a peak signal-to-noise ratio (S/N) of over 300 per 1 A resolution element in the rest frame. 
See Ref. 1247J for a comprehensive line list of 80 emission features. Similar composite 
spectra have been recently constructed from HST Faint Object Spectrograph and FUSE 
observations f209l 12901. The main advantage of composite spectra is that S/N is so high 
that identification of many faint features that remain invisible in individual spectra becomes 
possible \J1\. Most prominent features and overall continuum shape are easily appreciable; 
from Fig. ^one can realize the extent of Fell emission (shaded areas), the prominence of 
hydrogen Balmer (especially broad Ha, H/3, H7), of Lyman lines, and of the ClvA1549 
line, which makes them, apart from their physical importance, the best studied optical/UV 
lines of AGNs. Beyond line identification and some basic information a composite spec- 
trum can yield deceptive results. In principle it is legitimate to median together all spectra 
of a given survey if spectral properties scatter randomly with reasonable dispersion around 
the median. This is not the case for color selected samples of quasars. It is important to re- 
mark that attempts to model a composite spectrum like the ones from the SDSS and LBQS 
fTTl 12471 may be doomed to failure. Ionization conditions and gas kinematics are not the 
same for all type-1 AGNs 11391 . The use of average line ratios from composite survey 
spectra has led to an impasse in the attempt at reproducing the observations through sim- 
ple photoionization models (see e.g., Ref. 1 170 1, and references therein), an impasse that 
still hampers present-day efforts. Recent high quality data 12121 with nearly simultaneous 
FUSE, HST and optical observations emphasize the quasar diversity in terms of emission 
line and continuum spectrophotometric properties. Optical data for w 200 AGNs at z <, 0.8 
obtained with resolution A/ A A ~ 1000 show the impressive diversity in terms of emission 
line profiles 11401 . 

2.2 Basis of Eigenvector Analysis 

Eigenvector techniques are applied whenever many variables appear to be more or less 
loosely correlated without an intuitive indication of a dominant correlation or variable. A 
set of n objects may have m measured parameters like flux, FWHM, and equivalent widths 
of optical and UV emission lines, continuum shape, etc. We can think that each set of 
measurement is a vector in an m-dimensional space described by orthogonal vectors if, 
and define a matrix M of n vectors with m measurements. The Principal Component 
Analysis (PCA) searches for the best-fitting set of orthogonal axes to replace the original m 
axes in the space of measured parameters. The new axis set is sought by maximizing the 
sum of the squared projections onto each axis i.e., {Mlf)^{Mlf), where denotes the 
transposed matrix. If the set of measurements has been previously centered subtracting the 
variable average, then M^Af can be thought as a variance/covariance matrix. In spectral 
PCA no measurements are performed: the whole spectrum of n objects is divided into 
m small wavelength bins, and each input variable is the flux in the wavelength bin |213| 
|285l. We seek the maximum of if M^Mlf imposing the condition that the norm of if is 
unity through a Lagrange multiplier A, and setting the first derivative to 1 156,1 . With our 
formalism we have if"^ M'^ Mlf — X(lf^lf — 1), hence IM^Mlf — 2\lf = 0, and then 
M^^ Mlf = Xlf. This is an eigenvalue equation, which can be solved numerically. More 
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eigenvectors are sought through a similar procedure, with the additional constraints that the 
eigenvectors must be orthogonal to each other. Since the eigenvalues are a measurement 
of the sum of the squared projections of the original data vectors on the new normalized 
vectors, and since we assume to have set as a covariance matrix M'^ M, the eigenvalues are 
a measurements of the amount of variance carried in each new direction. For example, in a 
plane we can imagine a set of almost aligned points; their projections in the original axes 
may be nearly equal if they are not aligned preferentially with anyone of the two axes of 
the two original frame. In a physical context we may think of two variables that are highly 
correlated. We can maximize the projections along one axis by simply operating a rotation 
of the reference frame. A linear combination of the two original variables (or vectors) will 
now constitute the one vector that is needed to account for the variance of the data. A key 
aspect of the power of the PCA emerges from this simple example: a problem originally 
treated in two dimensions was inherently one-dimensional i.e., a PCA can restore a problem 
with a very large set of variables to its intrinsic dimensionality. 

2.3 The Original Eigenvectors by Boroson & Green 

The first important application of PCA to the interpretation of quasar spectra was done in 
the early 1990s on a sample of 87 objects of the PG survey at z < 0.5 l22ll . T. Boroson 
& R. Green measured the most prominent emission features in the H/3 spectral region, and 
found that most of the variance was related to two sets of correlations, the first being an 
anti-correlation between the prominence of [Olll]AA4959,5007 and optical Fell emission 
(Fellopt)- The "eigenvector 1" has been later found in a number of much larger samples 
ll2n l88l IT221 12291 12301 1285L and with spectral parameters describing phenomenologies as 
apparently distant as the X-ray continuum continuum shape and the radial velocity shift of 
high ionization emission lines. Even if many pieces of the eigenvector analysis had been 
hinted at by previous workers I18II138II263L the "eigenvector 1" provided, for the first time, 
a powerful systematic description of quasar spectral diversity. The second eigenvector in- 
volved the optical luminosity and the strength of the high ionization line (HIL) HellA4686. 
This second eigenvector had been also found by previous workers, and is basically the 
"Baldwin effect" ( ^Ot . 

2.4 The Optical Eigenvector 1 Plane 

It is remarkable that the spectral diversity of quasars can be reduced (if we exclude any 
luminosity dependence) to just two variables. The diversity and correlation of broad line 
AGNs are clearly shown by their distribution in the so-called "Eigenvector 1 (El) opti- 
cal plane," defined by the FWHM of the H/3 broad component, FWHM(H/3bc)> and by 
the equivalent width ratio between the Fe II blended emission at A4570 and the H/3 broad 
component, i?Fcii = W(Fell A4570)AV(H/3bc) 1121112211229 ^. The anti-correlation between 
FWHM(H/3bc) and -Rpcii provides perhaps a superior description of El than the original 
one, since it involves parameters related to the broad lines only ([Olll]AA4959,5007 are 
narrow lines affected by radio loudness; II230I '). The modest spread of data points drawn 
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Figure 2: The Optical Plane of the eigenvector 1 of type-1 AGNs; abscissa is the equivalent 
width ratio between Fell A4570 and H/3bc> ordinate is FWHM of H/3bc in km s~^. The 
plane has been subdivided in spectral type according to Ref. |230|. The dot-dashed line 
identifies Narrow Line Seyfert Is, defined by the condition FWHM(H/3bc)< 2000 km s^^. 
In the outlier regions, very few genuine sources (i.e., excluding those that are borderline 
because of errors) are found. Data are from Marziani et al. 11401 



from a sample of more than 200 type-1 AGNs allow the definition of typical spectral types 
covering a narrow range in FWHM(H/3bc) and i?Feii |230|. We set Ai?Fo = 0.5, and, 
for FWHM(H/3bc) ~ 4000 km s^\ we define spectral type Al, A2, A3 in order of in- 
creasing i?Feii (see Fig. |2l Fell emitters are considered strong if -Rpeii ~ 1; extreme if 
Rpeii ~ 1-5). Sources with FWHM(H/3bc) ~ 4000 km s^^ very rarely show iipc ~ 0.5. 
If Rfc ~ 0.5, we define a second sequence of increasing FWHM(H/3) (A1,B1,B1+ . . . ) 
with AFWHM(H/3)f« 4000 km s^^. The spectral types are easy to identify observationally, 
even from a visual inspection of the optical spectrum (Fig. |3ll and most of them may have 
a direct physical meaning (§ IZIl ©■ FWHM(H/3bc) ~ 4000 km s"^ represents a 
notable limit. It separates two groups of AGNs, Population A and B I229II whose different 
properties are most likely related to structural differences in the Broad Line Region (BLR; 
see 



2.5 Eigenvector 1 Measurements 

Quoting Peterson et al. fV79\, we can say that accurate line-width measurement depends 
critically on avoiding contaminating features, in particular the narrow components of the 
emission lines. Fig. |4] illustrates the typical analysis and measurement procedure applied to 
optical and UV data, in the H/3 (left panel) and ClvA1549 (right panel) spectral region. Two 
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Composite spectra for different bins 




Figure 3: Spectral diversity organized along the El sequence for UV (left panel) and op- 
tical (right panel) spectra, for the spectral types defined according to Sulentic et al. 12301 . 
Prominent blends of Fell emission are shaded. 

sources - one representative of Pop. A and the other of Pop. B - are shown. After continuum 
subtraction, Fell can be measured through a template Fell spectrum. Measurement of Fell 
A4570 (i.e.. Fell emission integrated over the wavelength range 4434^684 A) parameters 
like flux and width are accomplished by constructing an array of template spectra within 
reasonable limits of scaling and broadening factors I22lll38[ll40l . starting from a template 
(the I Zw 1 Fell spectrum). One then obtains the best scaling and broadening factor by 
subtracting the templates and by identifying the one template that minimizes the sum of the 
least-square residuals in the range 4450^600 A. In principle. Fell and continuum should be 
fit simultaneously in an automatized way i75ll . especially if Fell is very strong. A similar 
procedure l59l BHl ITITl fT38l |258l is appUed to the ClvA1549 spectrum using an Felluv 
emission template 11381 12581 1276.1 . While the procedure works remarkably well in the 
optical and in the UV around 1600 A, some caution is needed in the UV between 2000 A 
and 3000 A since Fell emission around MgllA2800 can be strong and different from source 
to source 11041 . 

We usually isolate the contribution of H/3nc> [Olll]AA4959,5007, and of broad 
HellA4686 before we apply a high-order spline fit to produce a model- independent "de- 
scription" of H/3bc which minimizes the effects of noise (the thick lines in the side panels 
of Fig. m 1 140]). Measurements of FWHM(H/3bc) are one of the basis of black hole mass 
estimates, so that it is important to provide a reliable recipe that leads to accurate results. 
Subtraction of a narrow-line component H/3nc is best done according to the following cri- 
teria |140|: (1) if a clear inflection is seen, the subtraction is trivial. This is the case of 
most sources (see the example of NGC 3783 in Fig. |4ll; (2) in the case where an infec- 
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tion is not seen, we subtracted a Gaussian profile under the condition that FWHM(H/3nc) 
«FWHM([Olll] AA4959,5007). It is important to stress that this last condition is not always 
applicable; (3) for several Pop. A sources (e.g., "Narrow Line Seyfert Is" (NLSyls) like 
Ton 28 shown in Fig. HJi with typically Lorentzian profiles and FWHM ^ 4000 km s~^ 
( ^2.8I) . the inflection is not observed. Such sources are "blue outliers" i.e., objects for 
which the recessional velocity measured on [Olll] AA4959,5007 is lower than that of H/3nc 
by more than 250 km s^^ |288|. Subtraction of a significant narrow component in many 
of these sources would imply H/3nc stronger than or comparable to [Olll]A5007, which is 
not consistent with other forbidden line ratios 1 158 1. Typically, H/3nc is 1/10 the strength 
of [Olll]A5007. In the outliers, any H/3nc would be appreciably displaced along the H/3bc 
profile. However, since W ([Olll]A5007) ^ 2.5 A, the H/3nc becomes too weak to be 
detected (very high resolution spectra are needed; it is detected for I Zw 1 1255 1 . and ref- 
erences therein). This argues against the subtraction of a strong H/3nc in blue outliers if 
high-resolution does not resolve a low-ionization H/3nc [255 |. If this recipe is applied, 
FWHM(H/3bc) measurements with an accuracy of ±10% at a la confidence level are pos- 
sible for the wide majority of "good" spectra. We customarily isolate a narrow-component 
contribution also to ClvA1549 before the CivA1549bc is fit. Since the presence of signifi- 
cant CivA1549nc has been debated, we discuss the issue in physical terms in ^4.7.11 

2.6 Eigenvector 1 Correlates 
2.6.1 Fell 

The emission spectrum of Fell in AGNs spans a wide wavelength range from the UV to the 
IR. It affects not only line intensity measurements of other ions but the continuum determi- 
nation as well. In case of the strong Felluv. the emission produces a pseudo-continuum that 
extends to wavelengths as short as 1000 A, and is especially strong in the near UV between 
2000 and 3000 A (Fig. d 1 138ll250l 12761 '). A convenient subdivision of Fell emission may 
be as follows: (1) Felluv, between 2000 and 3000 A [104 250); (2) near Felluv between 
3000 and 3500 A; (3) Fell A4570, whose values are probably the most frequently reported 
Fellopt measurements [22 140 1. Under the assumption that the I Zw 1 spectrum is a good 
representation of optical Fell emission, it is Fell A4570/Fellopt ~ 3.3 lil38J : (4) Fell A5250 
i.e., the emission blend on the red side of H/3 (see Figs. Qandl^ll. 

Fell emission is still rather poorly understood, even if several recent papers have made 
significant advances after a decade-long standstill (see Ref. [229] for a bibliography before 
late 1999). An extremely complex energy level structure of Fell makes it very difficult 
to obtain all the experimental transition probabilities and therefore calculate line intensi- 
ties (l6l 122 1 1 12521 . and references therein). Although Z » Zq does not seem anymore a 
requirement for explaining very large Felluv/ MgllA2800 ratios 12501 . there is still no gen- 
eral consensus on the excitation mechanism. There are two lines of thought: Fe+ is mainly 
produced by photoionization; coUisional excitation and florescence produce the observed 
strength and multiplet ratio. A recent paper 12491 accounts for the observed Felluv/Fellopt 
and Felluv/ MgllA2800 intensity ratios. Large electron density He (UV multiplets are fa- 



Accretion onto Supermassive Black Holes in Quasars 



133 




Figure 4: Line profile analysis for two sources in widely different positions along the El 
sequence (Ton 28, a NLSyl, and NGC 783), for H/3 (left panels) and ClvA1549 (right 
panels). The upper panel shows the continuum subtracted spectra, with the Fell contribution 
highlighted. The middle panel emphasizes the H/3bc profile as well as the HellA4686 line. 
Note that the [Olll]AA4959,5007 lines of Ton 28 are significantly blueshifted with respect 
to the rest frame; Ton 28 is a "blue outlier". 



vored because of the increase in the upper level population) and moderate micro-turbolence 
(broadening ~ 10 km s~^) can account for the combination of important diagnostics such as 
the Felluv/ MgllA2800, and Felluv/ Feiiopt intensity ratio 12501 . However, models based 
on photoionization are not void of difficulties: a recent work indicates that photoionized 
BLR gas cannot produce both the observed shape and observed equivalent width of the 
2200-2800 A Felluv bump unless there is considerable velocity structure corresponding to 
a microturbulent velocity parameter vturb ~ 100 km s^^, which strongly favor fluorescent 
processes [6|. A second line of thought ( I43II255I . and references therein) suggest that mod- 
elling the spectrum of the BLR requires a non radiative heating mechanism which increases 
the temperature in the excited HI region, thus providing the necessary additional excitation 
of the Fell lines I107II . The specificity of such a medium compared with a photoionized 
medium is its extremely low degree of ionization. The emission spectrum is made exclu- 
sively of low ionization lines (LILs) like Fell, Balmer lines and MgllA2800 for the electron 
density rie typical of the BLR. The photoionization models require very high ionizing pho- 
ton fluxes ~ 10^'' — 10^^ cm~^s^^ to explain very large Felluv/ MgiiA2800 ratios, while 
coUisionally heated media require a shielding mechanism to hide the emitting regions from 
the strong AGN continuum. On the basis of our current understanding of the line emitting 
region structure it is not clear where this shielded medium could be l43l l6ll. 

Even if the poor understanding of Fell processes hampers the use of Fell emission as a 
diagnostics, there are several important constraints related to El: 
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Figure 5: The left panel shows three examples with widely different FWHM(Fell A4570); 
it is shown that the I Zw 1 template, opportunely scaled and broadened, provides a good 
fit to all three sources (it provides a satisfactory reproduction to Fellopt of almost all 
sources II140I . The middle and right panel show the correlation between FWHM(H/3bc) 
and FWHM(Fell A4570), for FWHM(H/3bc)< 4000 km s^^ and FWHM(H/3bc)> 4000 
km s"^ respectively. 

• Fell A4570 emission is well reproduced, in a sample of 300 AGNs 11401 11431 by a 
scaled and broadened I Zw 1 Fellopt template. All sources did not show any strong 
deviation from the template within the limits imposed by S/N and resolution. Mea- 
surements of Fell width and strength are intrinsically difficult. Non optimal intrinsic 
width, spectral coverage, and S/N imply that only an upper limit to the equivalent 
width can be measured for faint Fell sources (see Fig. 

• Fig- El shows that the template method is able to reproduce Fell in cases of widely 
different width and strength. At zero order. Fig. |5l confirms that Fell and H/3bc are 
approximately of the same width, suggesting a common origin for these two LILs 
and confirming that -Rpeii is a meaningful parameter. There is however a difference: 
while for FWHM(H/3bc) ~ 4000 km s^^ we have FWHM(H/3Bc)~FWHM(Fell 
A4570) in a rigorous statistical sense, a sign test suggests that FWHM(H/3bc) ~ 
FWHM(Fell A4570) if FWHM(H/3bc) ~ 4000 km s^^l 142 1; 

• the Felluv/ MgllA2800 ratio seems to be El dependent, in the sense that sources with 
the narrowest FWHM(H/5bc) |140| show the largest values of Felluv/ MgllA2800 
I104in?l |. The Felluv/ MgllA2800 ratio increase along the El sequence may imply 
higher iron abundance but is also consistent with an increase in Uq, if density changes 
from Ue ^ 10^° cm"3 iqU ^m-s (250). Largest Felluv/ MgiiA2800(ss 
10) ratios require nc~ 10^^ cm^'^. Higher density is supported by several lines of 
evidence {^UM; EU- 
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2.6.2 Prominence of [Oiii] AA4959,5007 

The prominence of [Olll]AA4959,5007 is one of the most basic correlates along El, es- 
pecially if the parameter actually measured is the peak height over continuum 11211 l22l . 
Although not customary, this parameter yields a better correlation since profiles tend to 
broaden at low intensity 12881 . making W([Olll]AA4959,5007) a worse correlate. Emis- 
sion of [Olll]AA4959,5007 is not straightforward to interpret. If a sample covers low- 
z sources and spans a large z range, [Olll]AA4959,5007 measures are prone to aperture 
effects since the [Olll]AA4959,5007 emission morphology can be anisotropic |98|. Re- 
solved "narrow line regions" (NLRs i.e., where [Olll]AA4959,5007 is produced) reveal 
that line emission is not only anisotropic but also strongly influenced by radio jets, in 
nearby RQ Seyfert nuclei and in BLRG alike 1 15. .36. .70. .271 J . Not surprisingly, the 
EW of [Olll]AA4959,5007 is in RL than in RQ sources of the same optical spectral type 
llT25ll2n[T43l . In addition, [Olll]AA4959,5007 profiles can be strongly asymmetric with 
a predominance of blue-ward asymmetries (|271 1, and references therein). In the El con- 
text we find that large W([Olll]AA4959,5007) sources can have a blue-ward asymmetric 
wing; low W([Olll]AA4959,5007) sources (located in the lower left part of the optical El 
diagram; |288|) have emission entirely ascribable to a semi-broad, blueshifted feature that 
may resemble the blueward wing observed in stronger [Olll]AA4959,5007 emitters. 

2.7 Eigenvector 1 Extremes 

Narrow Line Seyfert Is NLSyls have been recognized as a distinct type of Seyfert nuclei 
since 1985 1 172 1. The defining criterion is that the width of the Balmer lines be less than 
2000 km s"^ The [Olll]A5007/H/3 ratio (but H/3 is not necessarily H/3nc!) is often smaller 
than 3 l46l 12151 . NLSyls usually (but not always) show rather strong Fellopt emission. 
Their soft X-ray spectra are very steep I248l l263ll and variable (e.g., [90"2411). The original 
definition is not void of difficulty, since: (1) a sharp discontinuity in many properties of 
type-1 sources is seen at FWHM(H/3bc)~ 4000 km s^^, and not at FWHM(H/3bc)~ 2000 
km s"^ {^R. In other words, in the range 2000 km s"^ < FWHM(H/?bc) ~ 4000 
km s~^ sources show properties consistent with the ones of NLSyl, although less extreme. 
(2) Sources with iircii ~ 0.5 (spectral type Al; |230|) may not be true NLSyls since they 
often show prominent [Olll]AA4959,5007, and their H/3bc profiles can be different from 
the typical NLSyl and bin A2 and A3 profiles; several sources are core-dominated RL 
sources and some belong to the second extremal type mentioned in ^1.21 With these caveats 
in mind, NLSyls seem to be drivers of the El correlations since they show and minimum 
FWHM(H/3bc) and may have extreme i?Feii in the "main sequence" of the optical El plane 
II229I . In the context of El they are interpreted as young or rejuvenated gas-rich objects with 
low black holes mass and high accreting rate l46ll88i rri9l[T44llT8()l . 

BAL QSOs High-z and low-z BAL QSOs share some of the properties of the extreme 
Pop. A sources in the El optical diagram: weaker CivA1549bc in emission l97l 12421 : 
stronger Felluv and AI111AI8O8 I261I . as well as stronger NvA1240. BAL QSOs occur 
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Figure 6: Inter-profile H/5 and ClvA1549 comparison for the prototypical BAL QSO PG 
1700+518. The profiles are shown here after continuum and Fell subtraction. Note that the 
ClvA1549 emission component is almost fully blueshifted with respect to the quasar rest 
frame. 



with a frequency of 15% with respect to all quasars 11961 : the fraction of radio-loud 
BAL quasars may be consistent with the fraction of radio-loud ordinary quasars f2391. The 
"Balnicity" Index |269| distribution rises with decreasing index values. This suggests that 
the fraction of quasars with intrinsic outflows may be significantly underestimated 1 196 1. 
Several high-z BAL QSOs have extremely strong Fell and very weak [Olll]AA4959,5007, 
extending the inverse relationship found for low-z QSOs and typical of Pop. A |287|. 
An important information that is still missing for most high-z BAL sources is a reliable 
estimate of the systemic redshift I117II . Without this information, the interpretation of the 
ClvA1549 absorption/emission profile is rather ambiguous. Identification of low-z BAL 
sources had to await space-based UV observatories 1231 12251 1243 i . and only a handful of 
sources is known to-date. Classical BAL sources with Balnicity Index » km s^^ at 
low and high z for which the rest frame can be reasonably set from narrow lines show an 
almost fully blueshifted emission component of ClvA1549 in addition to a fully blueshifted 
absorption (|225|; Fig. |6ll. Most low-z BAL AGNs with terminal wind velocities in the 
range « 20000-30000 km s^^ are either outliers or are located along the upper envelope 
of the El "main sequence" for Pop. A sources. Most revealing are the cases of Mkn 231 
and IRAS 07598-1-6508, which are basically the only 2 outliers in a sample of ~300 low-z 
AGNs. The sources show FWHM(H/3bc) significantly larger than FWHM(Fellopt). The 
difference is due to a strong blue-ward asymmetry visible in the H/3bc profile. The H/3bc 
profile can be interpreted as due to a narrower component with FWHM FWHM(Fellopt), 
and to a completely blueshifted feature, most probably associated to an high ionization 
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outflow emitting most of HILs like ClvA1549. If this effect is taken into account, and 
FWHM(Fellopt) instead of FWHM(H;3bc) is used, Mkn 231 and IRAS 07598+6508 move 
into the "main sequence" optical El diagram to become extreme (in terms of i?Feii) Pop. A 
sources. 

Radio Quiet & Radio Loud Quasars FWHM(H/3bc) in RL sources are at least twice 
as large as the RQ majority |229|. The average broad Fell A4570 emission line strength 
is also about half that for RQ sources 1221 II06I 138 1. Fig. shows that the Fanaroff- 
Riley II (i.e., lobe-dominated, LD) "parent population" of relatively un-boosted RL sources 
(median radio/optical flux ratio Rk ~ 500) shows the most restricted occupation. The 
Doppler boosted core-dominated (CD) RL sources (median Rk ~1000) lie towards smaller 
FWHM(H/3bc) and stronger Fell in El as expected if the lines arise in an accretion disk. 

2.8 Additional El Correlates 

H/?BC Profile Shape Median spectra of the H/3bc profile were extracted for each spectral 
type (Fig. ^ after Fellopt and continuum subtraction in a sample of ~ 200 AGNs [ 140| 
|230|. The median H/3bc profiles in spectral types A are almost symmetric, with a slight 
blueward asymmetry in bin A2 (see Fig. |2lfor the definition of the spectral types). The blue 
asymmetry becomes visually apparent in bin A3 which contains sources that are extreme 
in many ways. The median profiles of Pop. B sources are very redward asymmetric, with 
the strongest asymmetry in the (few) bin Bl sources. The H/3bc profile in NLSyl and 
Pop. A sources is well fit by a Lorentz function [ 254J . The best fit to H/?bc bin Bl and 
B1+ profiles is achieved through the sum of two Gaussians: (1) an unshifted Gaussian core 
(FWHM(H/?Bc) ~ 5000 km s"^) and (2) a broader redshifted Gaussian base with FWHM~ 
10000 km s^^ 129,1. We remark that a Double Gaussian fit to H/3bc in Pop- B sources is 
a formal result. The two components have a physical justification if, for the generality of 
AGNs, the broader one can be ascribed to gas that lies closest to the continuum source in 
an optically thin (to the Hi ionizing continuum) region with large covering factor /c — > 1 
Il29ll30l TO 12141 . 

[Olll]AA4959,5007 blueshifts The narrow lines of [Olll]AA4959,5007 usually provide 
a reUable measure of the systemic, or rest-frame, quasar velocity. However, several ob- 
servations indicate that the NLSyl galaxy prototype I Zw 1 shows an [Olll]AA4959,5007 
blueshift of -500 km s^^ relative to other rest-frame measures l24l 1138 1. with a sec- 
ond component at -1500 km s~^ (see references in Ref. 1.255 J ). I Zw 1 is not alone. 
It is possible to identify sources with [Olll]AA4959,5007 Avr ^ - 300 km s~^ with re- 
spect to H/3nc- Such large shifts are unlikely to be due to measurement errors [288 1. The 
two known largest blueshift after I Zw 1 reach Av^ ~ —1000 km s^^ |2|. Since these 
sources are apparently rather rare (but not unique! 189l|9l][l4l), and since they seem to 
lie out of a continuous shift distiibution, they have been dubbed "blue outliers". The dis- 
tribution of blue outliers in the El optical plane is different from that of the general AGN 
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Figure 7: Distribution of RL sources in the optical plane of the eigenvector 1 of type-1 
AGNs; abscissa is the equivalent width ratio between Fell A4570 and H/3bc> ordinate is 
FWHM of H/3bc in km s^^. The big arrow indicates displacement between the median 
Rpcii and FWHM(H/3bc) values from lobed and pure core RL sources. The circle identi- 
fies CD sources that genuinely have -Rpeii ~ 0.5 12341 : other sources at -Rpeii ~ 0.5 are 
strongly affected by measurement errors. 

population 12881 . The blue outliers occupy the lower right part of the diagram and are 
exclusively Pop. A/NLSyl nuclei [289|. Low W([Olll]AA4959,5007) may show only a 
blue-shifted component as if only the asymmetric part of the profile observed in stronger 
sources is being emitted. It is relatively straightforward to interpret the large blueshift and 
the profile of the [Olll] AA4959,5007 lines as the resuh of an outflow l2l l288l . and the low 
W([Olll]AA4959,5007) as a signature of a compact NLR I 255ii288l . 

High Ionization Lines The peak shift of the prototypical HIL ClvA1549 is an important 
El correlate in the sense that large shifts are shown only for FWHM(H/3bc) ~ 4000 
km s-i |4|[T38l|23l| at low z. W(ClvA1549) decreases along the El sequence from sources 
with large FWHM(H/3bc) and low i^pcii (typical W(CivA1549bc)~ 100-^200 A) to the 
most extreme NLSyls (typical W(CivA1549bc) ;S 50 A; I4. .138..2?)3ll23n \ 

Ionization Level/Electron Density A decrease in equivalent width of Civ A 1549 and 
H/3bc along with an increase in prominence of Fell A4570 in the El sequence toward 
NLSyl sources suggest a systematic decrease in ionization level from RQ Pop. B to 
Pop. A |139|. The behavior of W(ClvA1549) and W(H^bc), and of the diagnostic ra- 
tios Silll]A1892/ Clll]A1909 and AI111AI8O8/ Clll]A1909 can be understood in terms of 
a decrease in ionization parameter U (defined as the ratio between the ionizing photon 
flux and the electron density Ue) and an increase in rio from log [/ ~ — 1, log no~ 9.5, 
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to logU ~ -2 ^ -2.5, log nc~ 10.5 11 miTTllIlll. Notably, the diagnostic ra- 
tios indicate that log rig ~ 10.5 -^11 towards the NLSyl domain (see also Ref. L255ll 
for I Zw 1). 

X-ray Continuum The soft (energy ;$ 2-3 KeV) X-ray excess l45l |2721 is perhaps the 
most distinctive X-ray feature that differentiate sources along El f263"229|. The nature of 
the soft X-ray excess is still debated, although direct thermal emission from the accretion 
disk seems to be ruled out on the basis of the very high temperature derived for this optically 
thick component I186II . Even the newest X-ray data cannot constraint the current models 
I28l l85l. In the context of El several recent studies (including XMM-Newton observations; 
P186' '189 1) reveal systematic differences in the X-ray properties of RL and RQ quasars. 
XMM observations 1 1 86 1 confirm a previously reported anti-correlation between soft and 
hard photon indices and FWHM(H/3bc) I195[l263i . The soft X-ray spectral index correlates 
with the hard-X ray spectral shape, in the sense that NLSyl and Pop. A sources have a 
steeper spectral shape 1 129|. 

UV Continuum Inferences on UV/FUV continuum shape are complicated by the lack of 
data at low z before the advent of HST and FUSE, as well as by the possibility of strong 
extinction. A breakthrough occurred with the discovery in the spectrum of 3C273 of an 
optical-UV feature called the "Big Blue Bump", corresponding to an important fraction of 
the bolometric luminosity. Some soft X-ray selected AGNs are NLSyl galaxies with an 
enhanced big blue bump emission component relative to the underlying continuum mak- 
ing the optical continuum blue and the soft X-ray spectrum steep i90ll . On the other hand, 
several recent works indicate that NLSyls have a redder continuum than the rest of type- 
1 AGNs 152. .46. ,2661. After correcting with an extinction curve obtained from ~1000 
AGNs, it seems that radio-quiet AGNs have very similar intrinsic UV-to-optical continuum 
shape over several orders of magnitude in luminosity, and that that radio-loud and radio- 
quiet AGNs probably share the same underlying continuum shape f^3). A recent investi- 
gation of ultraviolet-to-optical Spectral Energy Distribution (SED) of 17 AGNs used quasi- 
simultaneous spectrophotometry spanning 900-9000 A (rest frame) with FUSE, HST and 
the 2.1-meter telescope at KPNO to study the big blue bump. Most objects exhibit a spec- 
tral break around 1100 A, with a slope slightly steeper for RL sources [^ l2T2ll2?7ll^0l . 
but otherwise there is no obvious trend which may be related to the source location in the 
El sequence (apart from extinction). 

FIR excess Several quasars show an SED with an extremely strong FIR excess l20l l26l 
I2021I225I . The strong excess can be due, at least in part, to the contribution of circumnuclear 
starburst, as in the case of Mkn 231 |57|, or to reprocessed radiation from a dusty torus 
encircling the accretion disk and the BLR. NLSyls and extra-strong Fell sources are known 
to be found often among sources detected by the Infra-Red Astronomical Satellite 11331 
I154L although it is as yet unclear how frequent is a FIR bump among NLSyls 12151 . 
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Bolometric Correction It is not easy to obtain Lboi for individual quasars because they 
emit significant power over a large part of the electromagnetic spectrum, most notably in 
the FUV and in the FIR. Bolometric correction factors are found to be Lboi~ 9 13AL;^ 
ergs s~^ for A « 5400 A, with possible deviations for single sources of amplitude ±50% 
I65ll212ll281l . A bolometric relationship Lboi= 9^WXLx ergs s~^ has been often adopted 
EH [Toll [1431 [262 . 281 1. Recent work based on good coverage of the UV, FUV and X but 
still poor sampling of the IR SED confirms that the assumption L = 13AL;^ ergs s~^ at 
5400 A is a fairly good approximation |212|. Again, scatter is pretty large (±50%). To 
test any dependence of the bolometric correction along El, we analyzed 44 sources that 
share accurate H/3bc measurements 1 140| and SEDs constructed from archival data |281 1. 
Comparison of Lboi values computed from the SED and from the bolometric correction 
showed no significant systematic difference: the average AlogLboi is ~ —0.06 for all 44 
sources, with a standard deviation a 0.24. Similar results hold if RQ Pop. A sources and 
RL AGNs are considered separately. It is worth noting that the scatter is due to a minority 
of badly behaving data points (f«20%). If they are removed, the standard deviation becomes 
a 0.1 in all cases, with systematic differences always ^ 0.05. Although big blue bump 
studies show that there may be a gross uniformity as far as the FUV is concerned, the 
difficulty to assess prominence and inherence of any FIR bump as well as the still uncertain 
shape in the EUV warrant caution before claims of accuracy better than ±50 % can be 
believed for the bolometric correction of individual objects. It is also important to stress 
that sources which are RL CD could add an appreciable scatter to an average SED, since 
their continuum may be strongly affected by relativistic beaming. RL CD sources should 
be considered on an individual basis and not included in averages. 

X-ray Variability A common phenomenon observed in AGNs is variability (for a review 
see e.g., Ref. |244|). The most vigorous variability observed is that of the X-rays [1571. 
Giant-amplitude X-ray variability by more than a factor of 15 has been found in NLSyl 
galaxies. The most extreme one is IRAS 13224-3809 which raised by a factor of about 57 
in just two days fT9ll . This giant variability is persistent f79l I286II . Examination of X- 
ray properties across the Seyfert population reveals an anti-correlation between variability 
amplitude and FWHM(H/3bc) 191 2411. 

Optical Variability Optical variability has been established as an identifying property 
of type-1 AGNs for more than three decades. Typically these objects show continuum 
variations by 1-2 magnitudes with timescales ranging from days to years P185','267 1. Broad 
emission lines have also been found to vary. Emission line variations lag the continuum 
variations with delays ranging from a few days to months ISTl I183II . The conclusion of 
an optical photometric study on six NLSyls is that, as a class, there is no evidence that 
NLSyls behave any differently than broader-lined Seyfert Is IllOI . However, an earlier 
study 1 268 1 conducted on a large sample (23 AGNs) found some evidence that AGNs with 
the strongest Fell emission lines were less likely to show strong optical variability. In the 
El context, an unpublished analysis of a survey of PG quasars l86l suggests that strong 
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optical variability Am 0.5 is not present among sources with FWHM(H/3bc) ~ 4000 
km s~^, and i?Feii ~ 0.5; it seems to be frequent for FWHM(H/3bc) ~ 4000 km s~^, and, 
interestingly enough, for sources with FWHM(H/3bc) ~ 4000 km s~^ and i^Feii ~ 0-5. 

Summing up, along the El sequence we see extreme X-ray variability in correspondence 
of the high i?Fcii> low FWHM(H/3bc) end; optical variability, on the other hand, seems to 
follow the opposite trend: to be modest or even undetected in NLSyls sources, and to be of 
relatively large amplitude for some sources with broader lines. 



Table 1: Main Trends Along the El Sequence 
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Radio 
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RQ and RL 
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ionization 


lower 
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3 Inferences on the Broad Line Region Structure 

Two observational strategies have been proved to be very powerful to investigate the struc- 
ture of the broad line region: reverberation mapping 1 1001 11761 and statistical inter-profile 
comparisons, especially between HILs and LILs l 3hl50lfT38ll . Two dimensional reverber- 
ation mapping (i.e., mapping done considering narrow radial velocity intervals in the line 
profile) requires special observational capabilities, and attempts to perform such measure- 
ments have yet to provide convincing results on BLR structure IIOOI . The correlations of 
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El already allow to organize data in a fashion that provides also appropriate input to rever- 
beration mapping studies. Table 1 summarizes the main trends along the El sequence from 
large FWHM(H/3bc) and low i?Feii, to the narrowest FWHM(H/5bc) and largest -Rpeii- 

3.1 A Different BLR Structure Separates Two AGN Populations 

An obvious question is whether there is a continuity in properties all across the sequence in 
the optical El plane. This does not seem to be the case. Three results stand out: 

1. a discontinuity cannot occur at FWHM(H/3bc)~ 2000 km s^^, the formal limit of 
NLSyls, since sources show continuity in properties at least up to FWHM(H/3bc) ^ 
4000 km s~^ L229J . For example, ultra-soft sources are found at FWHM(H/3bc)~ 
3000kms-i l27l : 

2. at FWHM(H/?Bc)~ 4000 km s^^ we observe a rather abrupt change in H/3bc profile 
shape; 

3. large CIVA1549 and [Olll]AA4959,5007blueshifts occur for FWHM(H/3bc) ~ 4000 
km s~^. Fig. ISlshows a correlation between FWHM(CivA1549bc) and ClvA1549 
peak shift that is valid for RQ Pop. A sources only; the same variables for Pop. B 
sources would show a scatter plot. 

These results have strong structural implications for the BLR. The simplest model 
able to explain H/3bc and ClvA1549 properties for sources with FWHM(H/3bc) ~ 4000 
km s~^ is a disk and wind system. An optically thick, geometrically thin disk emits most 
of H/3bc and all of Fell; optically thin gas with an outflow radial velocity component Vout 
emits ClvA1549 and other HILs ( I25lll55lll90ll : due to the absence of intervening absorp- 
tion we may expect that the material is highly ionized). The optically thick disk obscuring 
the receding part of the flow and leaves a net blueshift in the CivA1549bc line profile 
f4. 62 138|. Other lines, like Silll]A1892 and Clll]A1909 are most probably emitted at 
the base of the flow in the LIL emitting regions as also suggested by some reverberation 
mapping results 11761 . This model introduces a dependence on the aspect angle 6 that can 
be strongly different for HIL and LILs. 

Observations first supported this scheme in the context of a RQ/RL comparison l62l 
I138II . It seems now that a more meaningful comparison is between sources with 
FWHM(H/5bc) ~ 4000 km s^^ (note that this boundary may be somewhat luminosity- 
dependent, ^8.4. It . and the rest of type-1 AGNs |229| because sources with 
FWHM(H/3bc) ~ 4000 km s~^ can be both RL and RQ, and their broad line spectra 
can be indistinguishable. The evidence in favor of this model is not only statistical but 
also based on detailed inter-profile comparisons and monitoring of individual sources. For 
example, NLSyls IRAS 13224-3809 and IH 0707^95 are characterized by very blue 
continua; broad, strongly blueshifted HIL ClvA1549 and NvA1240; narrow, symmet- 
ric intermediate-ionization lines (including Clll]A1909, Silll]A1892, and AI111AI8O8) and 
LILs Uke MgllA2800 centered at the rest wavelength 1 131 1. In NGC 4051, the HellA4686 
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Figure 8: Correlation between CivA1549bc shift and FWHM for radio quiet, Population 
A (FWHM(H/3bc ~ 4000 km s^^) sources. A formal least-square fit is shown. 

line is almost five times broader than H/3 and is strongly blueward asymmetric. Variabil- 
ity and single-epoch data are consistent with the B aimer lines arising in a low-inclination 
(nearly face-on) disk-like configuration, and the high-ionization lines arising in an outflow- 
ing wind, of which the near side is preferentially observed I180II . Akn 564 is the most 
extensively monitored NLSyl galaxy in the UV |40]. Absence of response in the canonical 
HIL ClvA1549 line is consistent with matter-bounded emission. 

3.1.1 Extreme BAL QSOs as Extreme Pop. A Sources 

Assembling the previous results with the main constraint set by much earlier work (see, e. 
g., Ref. l64ll ') can lead to a simple geometrical and kinematical model which encompasses 
Pop. A sources as well as BAL quasars (see Fig. |9ll. The full blueshift of the emission com- 
ponent of ClvA1549 (as in Fig.|6ll requires that the opening angle of the flow is ^ 100°. 
Closer to the boundaries of the outflow and beyond the BLR proper extends the BAL re- 
gion in a conical corona of divergence angle 10°. Whenever a secondary absorption is 
present (Fig. |6l see also Ref. 11171 for an atlas of BAL profiles), we may observe an axial 
region covering all of the continuum-emitting region and part of the BLR as well |225|. 
The depth of the absorption implies fc ^ 0.5. It is important to stress that a shell of absorb- 
ing material which could provide an adequate fc is not viable in the geometrical context of 
our model because of the absorption relatively narrow width. Bent flow lines seen at large 
viewing angle or along the flow could give rise to NAL and BAL respectively \64\. They 
would require an opening angle 180°, which does not seem supported by the ClvA1549 
profile in our sources. While the cylindrical sheet may appear rather ad hoc, it may have 
a straightforward physical explanation if it is associated to the axial flows which are likely 
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Figure 9: Schematic view of the central engine of extreme Population A sources. The 
region co-axial with the disk is postulated because of the presence of a secondary, broad 
absorption in most classical, low-z BAL QSOs. It may be present only in sources with 
significant black hole spin. 



present if the black hole spin is significant (©. 



3.1.2 On the Structure of Population B Sources 

The observational structural constraints are by far more ambiguous for Pop. B sources 
where we do not find a kinematical decoupling between HILs and LILs |62 1. This is not to 
say that a disk and wind system is not applicable, but that the evidence in favor of a wind is 
less compelling, probably because the ratio between outflow velocity and rotational veloc- 
ity Wout/^rot is much lower than in Pop. A sources. It is possible to fit the HIL profiles of 
Pop. B source NGC 5548 with a wind model, with parameters consistent with reverberation 
mapping results l39l . Modest outflows (~ 1000 km s^^) are seen also in Pop. B sources 
like NGC 3783 I120II . A mild stratification can explain why higher ionization lines show 
broader profiles |142|: the FWHM of HellA4686, H/3bc> and Fell decreases in this order. 
A "split" BLR explains the same trend in Pop. A \4\ and why the HellA4686, Hel, and H/5 
lines respond to continuum change with increasing delay in Mkn 110 L114>.lllil . 



3.2 The BLR and the Accretion DiskAVind Paradigm 

Emission lines originating from a geometrically thin, optically thick disk |44l will show 
extremely small FWHM when observed face-on. Observed properties of CD RL sources 
suggest that we need a considerable velocity dispersion in the vertical direction to account 
for observed LIL widths (FWHM(H/5bc) ~ 3000 km s~^). A candidate for the line emitting 
region involves the outer, self-gravitating part of the disk l42lfT4ll . At some distance rsF 
the disk is expected to become gravitationally unstable and to dissolve into individual self- 
gravitating clouds or rings. For the face-on case, assuming orbital motion with Keplerian 
angular velocity (JIk). one can write: 
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FWHM(H/3bc) ^Av^ ly^KRsG, 

where At; is the vertical velocity dispersion, assumed to be proportional to the Keplerian 
velocity by a factor of u. A reasonable guess for u is about 0. 1-0.2. One can easily show that 
FWHM of 1000 km s'^ implies an rsF~ 5000 Rg, and 3000 km s"! implies rsF~500 
Rg. The Toomre stability criterion applied to a standard Shakura-Sunyaev disk [21 1] yields 
results dependent on the assumption of the dominant source of opacity, making the rh and 
Mbh dependence of rsF highly uncertain 1 14, 42|. It is however reasonable to conclude 
that rsF can be smaller in Pop. B sources by a factor ~ 10, and that this may leave a very 
"small" emitting surface for any standard optically-thick geometrically-thin disk. Part of the 
line profile may be produced in the fragmented disk if it is illuminated by a geometrically 
thick, hot inner region {Advection Dominated Accretion Flow, ADAF, or evaporated disk 

mi 

The clumpy structure expected for r ^ rsF may be resolvable through high-resolution 
spectroscopy. We note also that double-peaked profiles are expected to be resolved if the 
outer BLR radius is tblr ~ 10^ i?g. However, double-peaked line profiles from a disk can 
easily be turned into single-peaked profiles by the presence of a disk wind 11551 . Although 
Voat ^^rot> the outflow Velocity gradient is as large as the rotational velocity gradient. 
Since photons can escape much more easily along lines of sight with a small projected 
velocity, the resulting line profiles are single peaked with broad wings even though the 
emission comes from gas that is essentially on circular orbits 11551 . 



4 Black Hole Mass Determination 
4.1 The Virial Assumption 

One can estimate the mass of the supermassive black hole using FWHM(H/?bc) and a 
reverberation BLR "radius" tblr 11081 along with the assumption of virialized motions 
L123J . The virial mass is 



G ' 

where G is the gravitational constant. \f v = FWHM of a suitable line, / \/3/2 if the 
orbits of the BLR gas elements are randomly oriented. 

The first question is whether the virial assumption is consistent with the data. If the 
virial assumption is valid, the BLR dynamics should be dominated by the gravity of a 
central point mass. In this case the characteristic line broadening should correlate with the 
time lag for different lines. The virial relationship is only marginally consistent with the best 
time-delay data 1 179','1231. However, emission line profiles of Pop. A sources are relatively 
symmetric and smooth [ 140 |, and the optically thick part respond in a roughly symmetric 
fashion to continuum changes in 3C 390.3 12171 and NGC 5548 l2T8l . In the following, 
we will discuss the appropriateness of the virial assumption to derive Mbh considering the 
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correlations and trends along the El sequence of §|2j as well as the structural constrains on 
the BLR derived in§|3] 

4.2 Photoionization Mass 

The first estimates of Mbh were based on the rough similarity of AGN spectra, and on 
the consequent assumption of constant ionization parameter U or of constant product Uue 
(P62 | and references therein). A cumbersome evaluation of the ionizing luminosity is still 
needed with these methods. Even if these estimates are very rough, and the assumption of 
constant U is debatable (but see Ref. IISOI I. these studies provide a consistency check for 
photoionization models and reverberation mapping results 12621 . Further refinements may 
be possible when the behavior of U along the El sequence will be better understood and if 
Fell can be used to constrain ionizing photon flux and He 1.139. .250J . 

4.3 Mass Determination through Reverberation Mapping 

The distance of the BLR tblr can be deduced from reverberation data, most notably from 
H/?BC data. The cross correlation function between the continuum and the emission line 
light curve measures a time lag tj^ between continuum and line variations. The time lag, 
due to the light travel time across the broad line emitting region, yields an estimate of the 
'"BLR~ ctL- It is important to stress that the derivation of tblr follows from several as- 
sumptions: (1) the continuum emitting region is much smaller than the line emitting region; 
(2) observable and Hi ionizing continuum are related. This assumption seems to hold well 
since the monochromatic luminosity at UV selected wavelengths strongly correlates with 
''BLR L150 2561; (3) the light travel time across the BLR is a most important parameter, in 
the sense that it is shorter than the dynamical time (so that the BLR structure is not changing 
over the light travel time); (4) no dynamical effects of radiation are considered. Some pe- 
culiar structures in the line profiles may however derive from radiation pressure P137'|; (5) 
the line response is linear. The responsivity of the Balmer lines are generally anticorrelated 
with the incident photon flux. Thus, the responsivity vary with distance within the BLR 
for a fixed continuum luminosity and changes with time as the continuum source varies 
1116T. More technical problems involve unevenly sampled data, errors in flux calibration, 
normalization of spectra (different setup, aperture effects), and dilution by stellar contin- 
uum I176lll00ll . The coupled effects of a broad radial emissivity distribution, an unknown 
angular radiation pattern of line emission, and suboptimal sampling in the reverberation 
experiment can cause systematic errors as large as a factor of 3 or more in either direction 
Q23J. 

The line FWHM observed at a single epoch is not necessarily the best estimator of the 
gas velocity dispersion. The highest precision measure of the virial product is obtained by 
using the cross-correlation function centroid (as opposed to the cross-correlation function 
peak) for the time delay and the radial velocity dispersion of the variable part of the line 
1184111791 . The velocity dispersion for the variable part of the spectrum (which is due to 
portion of the BLR that is truly optically-thick) has now become available for 35 sources 
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11791 . With those data Peterson et al. P179| found that the random component in the error 
of reverberation-based Mbh measurements is typically around 30%. 

4.3.1 Extending Reverberation Mapping Results through the BLR Size-Luminosity 
Relationship 

It is now common to estimate Mbh by assuming that the BLR distance from the central 
continuum source is iblr oc (L5100)" and a = 0.6 - 0.7, as derived from the reverberation 
data by Kaspi et al. (| 108 1; ^4.3. It . We therefore can write the black hole virial mass as 
follows: MBHtx FWHM(H/?Bc)^(Ai5ioo)"> where L5100 is the specific luminosity at 5100 
A in units of ergs s~^A~^ 11431 . 

Any deviation from a = 0.6 0.7 has quantitative effects. The mass of luminous 
quasars at z ^ 0.4 has been computed by extrapolating this relationship, so that it is im- 
portant to stress that the relationship is based exclusively on quasars of z ^ 0.4 and that 
the high (and low) luminosity ranges of the correlation are poorly sampled. If we consider 
sources in the luminosity range 43.4 ^ logL/L© ^ 45 (i.e. where the sources in Ref. 
1 1081 show uniform luminosity sampling), the slope of the best fit is a = 0.8, and could 
easily be as high as a = 1 without increasing significantly the fit standard deviation. This 
case may be appropriate even for the PG quasar luminosity range 11391 . In dwarf active 
galaxies, a 0.5 seems appropriate i.e., dwarf active galaxies show larger BLRs than the 
values predicted by the tblr- I'Ly relation for more luminous AGNs [264 1 . One must re- 
main open to the possibility that 0.5 a 1, and that a might even be a function of L. 
Changing a implies an L-dependent change in mass estimates: the slope of the luminosity- 
to-mass relationship is affected as well as the location of points in the Lboi/AfBH vs. Mbh 
diagram in Fig. If a restriction is made to the most likely a range, a = 0.6 0.7, the 
effect on Mbh estimates is « 0.2 dex 11491 . 

Summing up all sources of uncertainty, individual Mbh obtained from single-epoch 
H/?BC observations and from the tblr- luminosity correlation seems to be affected by 
errors as large as a factor 2 — 3 at 1 o" confidence level L256..179J . 

4.4 Mbh Estimation through an Optically Thin VBLR 

Reverberation mapping-based Mbh determinations are probably affected by the 
non-negligible size of the BLR that is optically thick to the Lyman continuum, so that the 
derived tblr is not a very well defined quantity. The presence of high-velocity, optically 
thin line emission is likely rather common in AGNs (1*214' '2321 and references therein). 
Typical supporting evidence includes variability in the H/3bc line core, coupled with the 
absence of variability in the line wings, or strong response in HellA4686 without change 
in H/3bc- Observations of quasar PG 1416-129 revealed a large decline in its continuum 
luminosity over the past 10 years |232|. In response to the continuum change, the "classi- 
cal" broad component of H/3 almost completely disappeared (the flux decreased by a factor 
PS 10). A redshifted very broad component H/3vbc persisted after the demise of the broad 
component 12321 . In an optically thick medium the intensity of a recombination line is 
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Figure 10: Behavior of the intensity ratio HellA4686/H/3BC as a function of the ionization 
parameter U and electron density ng. 

governed by the luminosity of the ionizing continuum. If the medium is optically thin the 
intensity of the same recombination line is governed by the volume and density of the emit- 
ting gas and is not directly related to the luminosity of the ionizing continuum. The H/?vbc 
luminosity can be written as: 

L{Hf3yBc) = 47rrvBLR"-o^/i^'iaiAr, 

where ai and vi are the effective recombination coefficient and the line frequency for 
H/?, rvBLR is the distance from the central continuum source of a shell of density Ue and 
/c ~ 1. The unknown tvelr can be computed given the luminosity of H/3vbc after having 
inferred ric from the profile ratio in the line wings of HellA4686 and H/3 ( 11361 fTTSl : see 
Fig. ^|for the expected dependence of the HellA4686/H/3 intensity ratio). ryeLR is a well- 
defined quantity because Ar ^ r if the optical thickness to the Lyman continuum is less 
than unity at S> 10^° cm^^. In principle, if we could be sure that the wings of the Balmer 
lines are due to optically thin gas, the determination of rvsLR would be possible even from a 
single line profile observation of H/3bc» if the density and the covering factor are known. It 
id tempting consider the FWZI and the luminosity of the non- variable part of the line profile 
of H/3bc from reverberation mapping spectra. This approach could be attempted in the near 
future since excellent data are now becoming available from reverberation campaigns [ 179 1. 

4.4.1 Not Only Gravitational Redshift 

The amplitude of the redward H/?bc asymmetry observed in Pop. A sources seems to be 
mass dependent |143|. A factor f«6 increase in the redward displacement of the H/3bc 
centroid at 1/4 peak intensity [c(l/4)] means that the H/?bc asymmetry may be due to 
gravitational and transverse redshift |49|. A non-Doppler shift is due to a purely gravita- 
tional term « GMbh/^'blr- A second term is due to Doppler transverse shift, which is 




Log(electron density) 
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1/7, where 7 is the Lorentz factor. If gas motions are virial at tblrC^ Rg), the two terms 
yield a shift 

3 GM 
Az — ^ . 

2 C^rBLR 

The Pop. B sources considered by a recent study 11431 have been subdivided in nar- 
row ranges of Mbh and Lboi/MeH- In the range 3.5 ^ log Lboi/AfeH ~ 3.9, the resulting 
rvBLR values are 0.005 pc and 0.01 pc for log Af = 8 and logM = 9 respectively, if 
we take the c(l/4) value as a conservative estimate of the redshift. If we also model the 
VBLR as a gas shell (/c ~ 1) with optical depth to the Lyman continuum r ^ 1, CLOUDY 
f/Tl simulations show the shell emission falls far short in explaining the VBLR luminos- 
ity [1431. The difference between the expected and observed VBLR luminosity is largely 
a consequence of the small shell radius required to explain the large Avr in the c(l/4). 
We conclude that, even if c(l/4) is mass dependent, the c(l/4) shift amplitude cannot be 
explained by gravitational and transverse redshift alone 11381 . 

4.5 Mbh and Host Mass 

It seems that luminous (-24 > My > -28) quasars (both RL and RQ) are hosted in galaxies 
which are spheroidal or, at least, possess large bulges |73|. A correlation of nuclear black 
hole A/bh with stellar bulge velocity dispersion a^, is now well established in nearby galax- 
ies II72III66I . Supermassive black holes in galactic nuclei are though to be closely related, 
even in an evolutionary sense, to the bulge of the host galaxy. Reverberation-based Mbh 
estimates can be calibrated using the the correlation between Mbh and o"^, even if indirectly 
(i.e., one cannot use the same objects unless the luminosity drops so much that a* becomes 
measurable; I160I '). Seyfert galaxies (if we exclude NLSyls, I145II ') follow the same Mbh 
relation as nonactive galaxies, indicating that reverberation mapping measurements are con- 
sistent with those obtained using other methods 1 166|. Results based on the reverberating 
part of H/3bc U79 1 suggest that the systematic uncertainty in the Mbh is less than a factor 
of 3 I166II . The relationship between Mbh and a^, should be taken with care at low Mbh- 
NLSyls seem to be often host in dwarfish galaxies 11191 . Analysis of PG quasar obser- 
vations suggests a nonlinear relation between A/bh and bulge mass (Mbh~ -^^buig^*^^^) 
although a linear relation cannot be ruled out 1 1261. The mean MBn/Mbuige ratio may drop 
from 0.5% in bright (My ~ -22) ellipticals to 0.05% in low-luminosity (My ~ -18) 
bulges (see also Ref. 11451 ). 

Even more uncertain is the relationship between Mbh and bulge absolute magnitude. 
Seyfert galaxies are offset from nonactive galaxies but the deviation can be entirely under- 
stood as a difference in bulge luminosity, and not in Mbh; Seyfert galaxy hosts are brighter 
than normal galaxies for a given value of their velocity dispersion, perhaps as a result of 
younger stellar populations |160|. We indeed observe post-starburst quasars i.e., type-1 
AGNs that also display the strong Balmer jumps and high-order Balmer absorption lines 
characteristic of a very massive stellar population with ages ~ 100 Myr, even if they are a 
few percent of the quasar population 11751 . 
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4.5.1 Mass Determination from [Oiii]AA4959,5007 

In high-luminosity quasars, the relationship between Mbh and o"* can be studied compar- 
ing Mbh derived from the width of the broad H/3bc line and Mbr from the width of the 
narrow [Olll]AA4959,5007 lines used as a proxy to measure cr^. RQ AGNs seem to con- 
form to the established Mbh-c* relationship up to values of Mbh~ 10^° M©, with no 
discernible change in the relationship out to z w 3 [2201. There are two major difficulties 
here. Even if an A/bh - FWHM([Olll]A5007) correlation is present 1 161 1, scatter is large. 
FWHM([Olll]A5007) measures may not always provide a way to estimate Mbh- Mbh 
values from [Olll]AA4959,5007 can be considerably higher than values calculated using 
FWHM(H/3bc) II143I . This is not surprising because NLSyl-type AGNs apparently do not 
follow the same relationships as other type-1 AGNs 1 145 1. It is important to stress that low- 
W([Olll]AA4959,5007) sources can have FWHM([Olll]AA4959,5007) > FWHM(H/3), 
invalidating the virial assumption. Blueshifted [Olll]AA4959,5007 emission arises in out- 
flowing gas possibly associated to a disk wind |288|. The NLR in blue outliers may be 
very compact and its velocity field is not likely to be dynamically related to the host galaxy 
stellar bulge. This points to a hmiting W([Olll]A5007) below which FWHM([Olll]A5007) 
ceases to be a useful mass estimator. Only large W([Olll]A5007) RQ sources in Pop. B 
may have very extended NLR whose motions is dominated by gravity due to bulge stars. 

4.6 Masses of "Special" Sources 
4.6.1 "Double Peakers" 

A small fraction of AGNs exhibits exceptionally broad, double peaked LILs. The HctBC 
emission line profile is strikingly peculiar (see Ref. I229II for a few typical examples, or 
recent surveys such as the ones of Refs. |67 216|). AGNs with double peaked LILs remain 
relatively rare specimens, even if the SDSS has allowed to identify ~ 100 sources: w 4% 
(with many dubious cases) of all RQ and RL SDSS sources |216|. Although predominantly 
RQ, they are more likely to be found in RL sources, and account for about 20% of RL AGNs 
12161 l6Tl. Their relative rarity and their exceptionally broad lines prompted workers to 
search for a specificity either in terms of peculiar views or physics soon after the discovery 
of the prototypical source Arp 102B | 235|. 

Prototype "double-peakers" Arp 102B, 3C 390.3, and 3C 332 have been now moni- 
tored form more than 20 years I132II . A common property of double-peaked lines is slow, 
systematic variability of the profile shape on a timescale of years i.e., the timescale of dy- 
namical changes in accretion disks 167. 163 1 . Lower limits on the plausible orbital periods 
from the absence of peak radial velocity changes would require supermassive binary black 
holes |82| with total masses in excess of lO^'^ Mq. Such large Mbh values are difficult 
to reconcile with a maximum expected limit Mbh~ 3 • 10^ Mq (glU f66l), and with the 
known Mbh - host bulge mass relationship. Other recent results from a five year moni- 
toring of H/3 further support the dismissal of the binary black hole hypothesis for 3C390.3 
on the basis of the masses required 1217 .1. Hot spots, spiral waves and elliptical accretion 
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disks have been recently favored 11321122 41. If the hot spot lies within the accretion disk, 
it is possible to estimate Mbu and also the BLR physical size. The determination of tblr 
in physical units removes the degeneracy due to the rotational velocity field (i.e., the veloc- 
ity scales as (Mbh/ r)~^'^; disk model profiles yield a distance normalized to Rg). Mbh 
(2.2 lb 0.7) • 10^ Mq was inferred from the period and the distance of the hot spot in 
Arp 102B (4.8 • 10"'^ pc; [163]). It interesting to note that, within the framework of the 
model, Mbh/v 3 • 10^^ g cm~^ is still below the formal requirement for a black hole 
but the implied density p >, 10^^ Mopc^'^ is a stringent limit. The Mbh in double peakers 
can be also estimated using the empirical relationship between tblr and optical continuum 
luminosity 1 108 1. Mbh and m computations have been carried out for 135 objects from the 
SDSS |283| and from a survey of RL broad emission line AGNs l67ll . Mbh values range 
from 3 -lO"^ Mq to 5 -10^ Mq, and m is between 0.001 and 0.1. Double peakers are found 
up to L ~ 10^^ ergs s~^. 

4.6.2 Black Hole Binaries 

Several observational properties of extragalactic radio jets, such as bending, misaligment, 
and wiggling (often associated with knots superluminally moving along different-scale 
curved trajectory) have been interpreted in terms of helical structures of the jets. This 
structure is likely caused by the precession of the jet in a binary black hole system 1 35] 
|32li2l[lll]|260|. In the case of OJ 287 the bending of the VLBI jet was reported in Ref. 
II259I . A very small change in the orientation of the jet is needed in order to change the 
Doppler boosting dramatically, thereby producing long-term periodic brightness modula- 
tions. OJ287 was first recognized as a candidate binary black hole system by noticing the 
regularly spaced outburst pattern in its historical optical light curve 1 223 J . Photographic 
information on the brightness of this blazar (thought to be a variable star) extends for about 
100 years. Even though the observations were scanty in the beginning, there still was a 
convincing pattern which led to the prediction of the next major outburst in OJ287 in the 
fall of 1994 . For the first time a predicted cyclic phenomenon was indeed observed in an 
extragalactic object I222II . Sillanpaa et al. 12231 modelled the periodic outbursts using a 
binary system consisting of a pair of supermassive black holes with an orbital period of 8.9 
yr (in the rest frame of OJ 287). The light variations in this model are related to tidally 
induced mass flows from the accretion disks into the black holes 1 109,1. A variant of this 
model explains and predicts other features of the observed light curve 11281 . The variant 
allows for relativistic precession in the binary black hole system, and associates the major 
flares with the times at which the secondary black hole crosses the plane of the disk of the 
primary. The planes of the disks are perpendicular to each other. The monitoring campaign 
organized to observe OJ 287 in 1994 gave such a detailed light curve that a unique deter- 
mination of the orbital parameters was possible. The model gave masses of 10^ Mq and 
1.7 • 10^ Mq for the secondary and primary black hole respectively, and also predicted an 
eclipse of the secondary black hole disk by the disk of the primary in 1989. The eclipse was 
indeed observed 1246111941 in the optical but not in the radio, as expected. The observation 
of the next predicted flares for 2006 or 2007 (the date depends on the model) will allow to 
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Figure 11: Examples of sources with unambiguous CivA1549nc emission. The left panel 
shows two FWHM(ClvA1549) measurements: one for the whole ClvA1549 profile, and the 
second one for the CivA1549bc (thick solid line) only. Neglecting CivA1549nc results 
in large FWHM(ClvA1549), and hence Mbh measurement errors. The middle and the 
right panel show two radio loud sources with very high W([Olll]AA4959,5007) and widely 
different W(CivA1549nc)- 

refine models and make a direct measurement of the orbital energy loss in the system. 

4.7 At Low and High z: Mbh Determination from CivA1549 
4.7.1 Existence of a CivA1549nc Contribution 

The very existence of significant ClvA1549 narrow emission (CivA1549nc) has been a 
contentious issue. There is plenty of evidence that at least some AGNs with strong nar- 
row lines have a strong CivA1549nc- Perhaps two of the most striking examples are PG 
0026+12 (RQ; a borderhne NLSyl with strong NLR emission) and Pictor A (RL; Fig. 
fTTT i. These are two cases in which an obvious inflection between broad and narrow com- 
ponent is seen, and in which a CivA1549nc can be easily isolated with FWHM ^ 1000 
km s^^. The visual impression (already suggested in Ref. fTSSl) of a connection between 
[Olll]AA4959,5007 and the ClvA1549 narrow core is confirmed by a loose correlation be- 
tween W([Olll]AA4959,5007) and W(CivA1549nc) which we found considering only the 
cases of least ambiguous inflection. The data dispersion is probably intrinsic: at the two 
CivA1549nc extremes of large [Olll]AA4959,5007 we find Pictor A (W(CivA1549nc) ~ 
30 ±6 A) and 3C 390.3 (W(CivA1549nc) 5 ±1 A) with W([Oiii]A5007)~ 80 A (see 
Fig. fm . We do not expect a strong correlation physically, since emission of ClvA1549 is 
favored in the innermost regions of the NLR, as shown below. 

Even if RL sources often show prominent ClvA1549 cores of width ~ 2000 km s~^, 
such emission has been ascribed to CivA1549bc 1277 | and to an intermediate line region 
l30l . It is relatively naive to beheve that [Olll]AA4959,5007 and CivA1549nc should have 
the same width. Collisionally excited lines such as ClvA1549 show a band of efficient 
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reprocessing running at constant U along a diagonal ridge in the plane photon flux versus 
Tig. Emission is locally optimized along this ridge which corresponds to log U ~ —1.5 for 
ClvA1549 Q- Provided that a smooth density gradient exists (i.e., that there is a significant 
amount of gas at ne ~10^ cm~^), we expect ClvA1549 emission to be strong at larger rie 
(ClvA1549 emissivity is cx n1 while [Olll]AA4959,5007 becomes coUisionally quenched). 
If velocity dispersion increases with decreasing distance from the central continuum source 
it is somewhat natural to expect FWHM(CivA1549nc) > FWHM([Olll]AA4959,5007). 
Detailed calculations based on these simple considerations indeed suggest that 1000 km s~^ 
< FWHM(CivA1549nc) ~ 2000 km s"^ for a variety of possible density-gradient 
laws 1 170 226). The same models correctly predict that FWHM(H/3nc)~ 500 km s'^ > 
FWHM([Oiii]AA4959,5007) < FWHM(CivA1549nc)- To ascribe CivA1549 core emis- 
sion to the BLR (where rie ^ lO^'^cm"'^) seems rather arbitrary. Since CivA1549nc is 
prominent in Pop. B sources but can be absent in Pop A., the wind that gives rise to the 
blue-shifted CivA1549bc component may serve as a filter to the FUV ionizing radiation 
from the center that would otherwise reach the outer self-gravitating parts of the disk, where 
CivA1549nc may originate 11301 . 



4.7.2 Is the Width of CivA1549bc a Rehable Mbh Estimator? 



Estimating Mbh from the width of ClvA1549 is cumbersome. Failure to account for 
CivA1549nc has dramatic consequences on Mbh and Eddington ratio estimates (^. Dy- 
namically, the HIL emitting gas is not in virial equilibrium at least in a significant fraction 
of quasars (^ Fig. ISj. Assuming virial equilibrium for Pop. B (as suggested by the absence 
of systematic line shifts and by the tentative similarity with the H/3bc profile) still leaves the 
problem of properly assessing CivA1549nc- Recent studies of the ClvA1549 profile using 
HST archival spectra (|12, 122^ 11211 1256 , 266 1) have subtracted little or no H/3nc from 
the ClvA1549 line profiles. A comparison of line widths and shift measures for H/3 and 
ClvA1549 1 121 1266 1 shows that there are significant and systematic differences. An appar- 
ent dichotomy occurs at FWHM(H/3bc)~ 4000 km s^^ if no CivA1549nc is considered 
|FT2|: below 4000 km s^^, the ClvA1549 line is broader than H/3bc> but the reverse seems 
to hold when FWHM(H/3bc) ~ 4000 km s"^ where we believe that the NLR becomes 
important (Pop. B; |4|). This result is not necessarily against the view that ClvA1549 gen- 
erally originates closer to the center than H/3 since CivA1549nc was not subtracted. In 
the case of PG 0026-1-129, failure to subtract the NLR emission yields spurious broad line 
parameter measures (e.g. FWHM ~ 8000 km s^^ instead of 1860 km s^^; see Fig. fTTT i. 
We conclude that: (a) Mbh computations from ClvA1549 width are wrong in the case the 
ClvA1549 profile shows large blueshifts i.e., for Pop. A sources; (b) failure to properly 
correct for CivA1549nc yields a large error in the mass estimate |4|. The Mbh error is 
especially large for Pop. B sources; however, it must be remarked that CivA1549nc can be 
strong also in Pop. A RQ sources, and that a case-by-case analysis should be made. 
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4.8 High-^ Masses: Are They Really That Big? 

If the empirical relationship between tblr and the source luminosity is used to obtain 
Mbh at high z, the largest A^bh are ^ lO^'^ M©. Such values of Mbh suggest bulge 
mass ^ 10^^ M0 and cr^ ^ 700 km s^^ which are not observed at low z. Black holes 
with Mbh ~ 3 • 10^ M© should reside almost exclusively in high-redshift quasars 11691 . 
This implication is in contrast with the expectation that Mbh can only grow with cosmic 
time on timescales of the Universe present age, and inconsistent with several suggested 
scenarios of black hole and galaxy formation. A solution of this dilemma may reside in the 
improper use of ClvA1549 as a proxy of H/3 to estimate Mbh- A very good correlation has 
been recently found between tblr and the specific continuum luminosity at 3000 A 1 150|. 
In a sample of objects with broad-line radii determined from reverberation mapping the 
FWHM of MgllA2800 and H/3 are consistent with an exact one-to-one relation, as expected 
if both H/3bc and MgllA2800 are predominantly emitted at the same distance from the 
central ionizing source within a factor of 2.5 ri49l 11501 . The MgllA2800 line seems to 
be systematically narrower than H/3bc probably because MgllA2800 (like Fell and the 
reverberating part of H/3bc) is emitted in the most optically thick portion of the BLR. 
The FWHM(MgllA2800) can then be used to estimate Mbh of quasars 0.25 < z < 2.5 
via optical spectroscopy alone, and it is in principle even preferable to FWHM(H/3bc)- 
There are no strong theoretical reasons or empirical evidence to doubt the assumption that 
the MgllA2800 line broadening is virial, at least no more than the ones known for H/3bc- 
FWHM(MgllA2800) typically indicates a factor of 5 times lower Mbh than ClvA1549 
11601 . It is interesting to point out that a mass of 3 -10^ for the central black hole in the 
z = 6.41 quasar SDSS J114816. 64-1-525150 (one of the most distant quasars known) has 
been estimated applying the Mbh technique appropriate for a detected MgllA2800 FWHM 
of 6000 km s-i ||274l. This very high luminosity quasar does not show extremely broad 
lines and does not require super-Eddington luminosity. 

Measurements of H/?bc at high z can be achieved only thi^ough IR spectrometers. 
Instrumental capabilities lag behind the ones of optical spectrometers, and lack of resolution 
as well as poor SjN can lead to huge overestimates of FWHM(H/3bc) 11461 . although 
excellent data for a few tens of sources are becoming available 12331 . 

4.9 Mass Estimates Along the El Sequence 

The previous results suggest a number of general remarks. 

• The most accurate A/bh determinations are probably obtained considering the re- 
verberating part of emission lines 11791 . If Mbh estimates from the reverberating 
component of H/3bc are compared to ClvA1549 Mbh estimates, offsets of 1.0 dex 
or larger from a perfect one-to-one relationship are possible 1 256 1 . Mbh derived from 
CivA1549bc has a probability of just 90% to be within a factor 10 from the best re- 
verberation Mbh estimates. For the H/3 single-epoch estimates, the probability of 
getting A^bh accurate to within an order of magnitude is 95 %, which suggests 
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that systematic sources of uncertainties and ambiguities operate in the estimate of 
Mbh- 

• If an innermost region of optically thin (to the Lyman continuum) gas can be identi- 
fied from the profile of H/3, an estimate of the line emitting region size is relatively 
straightforward ( ^4.4t . However, the frequent occurrence of asymmetries partly in- 
validates the virial assumption for the innermost line emitting regions in Pop. B 
sources. 

• Single epoch, single line Mbu determinations based on the FWHM are valid in a 
statistical sense if the line is a LIL, like H/3 or MgllA2800. An estimation of the 
aspect angle ^ is a first-hand necessity, as farther outlined below. 

• It is improper to use ClvA1549 for Mbh calculations is the line is blueshifted or 
strongly asymmetric, as in Pop. A sources. This seems to be the case of most high-2; 
objects fl98|. TheMgiiA2800 line may provide the best estimator at high z. 

• Some special sources, like double peakers and black hole binaries may provide accu- 
rate Mbh. although tblr and A/bh determinations are still somewhat model depen- 
dent. Other sources - like BAL QSOs - may be favored if the orientation angle can 
be constrained. 

• [Olll]AA4959,5007 may provide a suitable Mbh estimator only in RQ sources, pre- 
dominantly for spectral types Al and Pop. B, if the recessional velocity from 
[Olll]AA4959,5007 agrees with the one deduced from other narrow lines which are 
of low ionization. 

4.9.1 Orientation Effects 

As reviewed in the previous sections, there is growing evidence that the LILs are predom- 
inantly emitted in a flattened system. LIL emission may be dominated by an extended, 
optically thick disk in Pop. A, while a fragmented disk may be the main line emitter in Pop. 
B sources. Inclination relative to a flattened system causes a systematic underestimate of 
the central mass. This effect is expected to be present in both RQ and RL sources. Consid- 
ering the vertical velocity dispersion for the line emitting gas at ^ ^ 0° FWHM(0 = 0), a 
suitable parameterization could be FWHM{e) ^ FWHM {6 = 0) + AFWHM ■ sin 6*, 
where FWHM{6 = 0) 1000 km s'^ and 3000 km s'^ for RQ and RL sources re- 
spectively. AFWHM can be estimated from the FWHM range from CD to LD sources 
( fl39. 234 1: a similar approach has been followed in Ref. 1 1481 ). Considering that 9 may be 
in the range few° 6 ^ 45°, and distributed randomly, we have an average < 6 >^ 30°. 
This means that ignoring orientation effects leads to a systematic underestimate of Mbh 
by A log Mbr^ 0.6. The underestimate may be a factor >, 10 if the source is observed 
almost pole-on. Even if pole-on sources are the rarest one in a randomly-oriented sample, 
errors of this amplitude may dramatically influence inferences on Mbh and blur correla- 
tions involving Lboi/MBH (which varies only by 2 - 3 orders of magnitude in AGNs; It 
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is also important to consider that optical luminosity may not be appropriate in the tblr - lu- 
minosity relationship for CD RL AGNs (9 0°) because jets may significantly contribute 
to the optical continuum. A relativistically boosted continuum leads to an overestimation 
of Mbh- In such cases, it may be better to consider an empirical relation between tblr and 
the H/3bc luminosity t284ll . 

5 Black Hole Spin & Observer's Orientation 

5.1 Where Are Spinning Black Holes? 

The idea of a non rotating black hole is rather distressing because every known massive 
object in the universe (from asteroids to neutron stars) does rotate. Conservation of angular 
momentum by the accreting matter should lead quickly to maximally rotating black holes 
|23 8 1. While Mbh is destined to increase on timescales comparable to the present age of the 
Universe, a black hole angular momentum can reverse and decrease through a diversity of 
mechanisms fSO" 101]: the collapse of massive gas accumulations, gas accretion, capture of 
stellar mass bodies, and successive mergers with other massive holes. A significant black 
hole specific angular momentum a is required for driving relativistic, radio-emitting jets 
(T7i although a maximally rotating black hole (a/M ^ 0.998G/c) is not necessary. In the 
El context, suggestions of a 7^ come from time scale variability arguments, especially 
from soft-X ray observations of NLSyl sources. A variation in luminosity cannot occur 
on an arbitrarily short time scale, a minimum physical limit being set by the light crossing 
time of the emitting region. The emitting matter should have some opacity r « nar ~ 1, 
where a is the Thomson scattering cross-section. The time needed for radiation to diffuse 
out of a region of size r is therefore At ^ {1 + T)r/c \69\. The change in luminosity AL 
following an accretion event with accretion rate AM/ At is AL = rjAM / Atc^ . Writing 
AM = ^TrR^nmp where n = r/aR and nip is the proton mass, and considering the travel 
time, we have: 

AL 4 nip 4^ o in4i -1 
At 6 a 

If changes of very large amplitude occur in such a short timescale that the relation above 
is violated for 7? w 0.1, claims have been laid that a rotating black hole is needed. This has 
been the case for a handful of NLSyl galaxies EH- In the view that some NLSyls and at 
least some other Pop. A sources are young/rejuvenated quasars II1441 12291 they may 
have experienced one of few accretion events leading to a consistent increase of the black 
hole angular momentum. 

Two-fluid models for relativistic outflows include a fast, relativistic beam surrounded 
by a slower, possibly thermal, outflow, with a mixing layer forming between the beam 
and the jet (e.g., [1341). The mixing layer could produce the secondary absorption seen in 
several BAL sources (Fig. |6ll. It is interesting to stress that the observed i?Fcii for the BAL 
QSOs Mkn 231 and IRAS 07598-1-6508 would require highly super-Eddington accretion 
unless the accreting object is a rotating black hole. Another aspect favoring a maximally 
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rotating black hole is the remarkable behavior of the Fe Ka emission line in MCG -06-30- 
15 12271 . This source is pretty unremarkable optically; it is a Pop. A source with modest 
Fellopt emission. More circumstantial evidence favoring rotating black holes is provided by 
the mass density of supermassive black holes in the present-day universe, which agrees with 
the quasai" luminosity evolution and total power output if the efficiency has been r] ^ 0.1, 
larger than that of a non-rotating black hole I 94J . 

It is then intriguing that, at present, neither black hole spin nor morphology seem to be 
sufficient conditions to explain radio loudness f63], disfavoring the idea that the RL and 
RQ dichotomy may be due to systematically different black hole spins [280 1. The physics 
of radio emission in the inner regions of all quasars (RL and RQ) may be essentially the 
same, involving a compact, partially opaque core together with a beamed jet lITOl . Other 
factors may be at play: a strong electromagnetic field, coupled with differential rotation, 
could serve to convert rotational energy into kinetic energy of outflows [34|. In the context 
of optical/UV observations such effects are at present observationally unconfirmed, also 
because radio loudness does not seem to have appreciable effects on the observed broad line 
spectrum 111431. In principle, we could expect that, if a viscous accretion disk is inclined 
in relation to the equatorial plane of a rotating black hole, the differential precession will 
produce warps in the disk. The combined action of the Lense-Thirring effect and the internal 
viscosity of the accretion disk forces the alignment between the angular momenta of the 
rotating black hole and the accretion disk at r ~ 1000 Rg ~ tblr- Some double peaked 
emission line profiles may be produced by warped disks. The agreement between some 
peculiar Hobc profiles and a toy model is very good L33. .3J. but many other models may 
also be applicable without invoking a warped disk 12241 . 

5.2 Can we Measure an Orientation Angle? 

In radio-loud quasars, where a rough orientation indicator is available from the radio core- 
to-lobe ratio, there is good evidence that FWHM(H/3bc) shows an orientation dependence 
consistent with a flattened (or even disk-like) geometry [275 . 291. There is a domain space 
separation between Fanaroff-Riley II (i.e., LD) and CD sources in El. The diagonal ar- 
row in figure Q indicates the average change in El position due to change in orientation 
from LD to CD sources. As stressed earlier, Mbh can be underestimated if no correction 
is applied to observed LIL widths 1 167 1. An analysis of the connection between Mbh and 
radio luminosity in radio-selected flat-spectrum quasars shows that values of Mbh are not 
systematically lower than those of luminous optically selected RL quasars, if a proper cor- 
rection for the effects of inclination is applied 1 105 1. It is also unlikely that there are sources 
radiating at Lboi/-^^Edd> 1 at z ^ 0.8 (see Fig. El 1 143 1) if a reasonable orientation cor- 
rection is applied to line widths. While a correction for the effect of varying 9 is feasible in 
a statistical sense, measuring ^ on a source-by-source basis has proved to be very difficult. 
Superluminal radio sources are an exception. Radio observations allow an estimate of the 
expected synchrotron self Compton flux relative to the observed X-ray flux. Lorentz factor 
and 6 can then be independently estimated ll84l . Figure[^shows 6 vs. FWHM(H/3bc) for 
II superluminal sources, confirming a factor « 3 effect on the H/3bc width at extreme 6 
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values. 




The observational properties of some, extreme NLSyls may well be due to a pole-on 
orientation angle. A NLSyl nucleus has been found in the extensively studied eruptive 
BL Lac, 0846+5 IWl, out of a large sample of NLSyl from the SDSS [292 1. The SDSS 
allowed the identification of another source with typical NLSyl properties (strong Fellopt, 
undetectable [Olll]AA4959,5007) that is definitely RL (Rk ^ 1000). The inverted radio 
spectrum and the very high brightness temperature derived from variation of the radio flux 
suggested the presence of a relativistic jet beaming toward the observer |291 1. It is unclear 
whether this finding can be extended to (at least) some RQ NLSyls. Some NLSyl might be 
considered the radio quiet equivalent of BL Lacs if they give rise to radio-silent relativistic 
outflow. Even if evidence is accumulating about the existence of a relativistic jet in RQ 
sources (e.g., 1 10 1), RQ sources (especially if of Pop. A) are not optically violently variable. 
Although some NLSyl sources (those with largest [Olll]AA4959,5007 and CivA1549bc 
blueshifts, as well as narrowest H/?bc and most strongly variable soft X-ray excess) are 
likely to be oriented almost pole-on, it is still not clear whether we are observing a beaming 
effect. 

Orientation estimates for RQ sources may come from the correlation between Fsoft 
and soft-X ray variability amplitude, although the poor understanding of the soft X-ray 
excess prevents the definition of an orientation indicator. Other attempts at estimating 9 are 
rather model-dependent. A sample of several thousand quasars from the SDSS confirms that 
HILs such ClvA1549 are significantly blueshifted with respect to LILs such as MgllA2800 
11981 . Among the SDSS sources, ClvA1549 emission-line peaks have a range of shifts 
from a redshift of 500 km s~^ to blueshifts well in excess of 2000 km s^^ compared to 
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Figure 13: Mass Luminosity relationship for a sample of k,210 AGNs [1431. Triangles 
and circles indicate radio-quiet and radio-loud sources, respectively. Open symbols indi- 
cate sources with FWHM(H/3bc) ~ 4000 km s~^ (Pop. B); filled symbols sources with 
FWHM(H/3bc) ~ 4000 km s^^ (Pop. A). The bars in the lower right corner indicate 
typical errors. 



MgllA2800. The anticorrelation between the shift of the ClvA1549 emission-line peak and 
W(ClvA1549) (Fig. ^ is confirmed by the SDSS study. Composite quasar spectra as a 
function of ClvA1549 shift suggest that the apparent shift of the ClvA1549 emission-line 
peak is not a shift so much as it is a lack of flux in the red wing for the composite with the 
largest apparent shift 1 198 1. An outflowing wind and an optically thick disk to act as a screen 
provide a tempting explanation. However, it is unlikely that purely geometrical effects are 
at play in determining the width and shifts of H/3bc and CivA1549bc alike 11391 . even if 
a simple toy model shows that ClvA1549 shift amplitude is orientation dependent with the 
largest shifts being produced by ^ ^ 0° |288|. W(ClvA1549) is strongly dependent on 
^boi/MBH (EIIH; see ^SATT i. as are the properties of outflowing winds |155 190]. To 
estimate on a source by source basis, we first need to understand how W(ClvA1549) and 
ClvA1549 shifts are affected by -Lboi/^Edd and 6. 
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We mention here an interesting attempt that is, again, strongly dependent on specific 
assumptions. Using the variable fraction of the H/3bc> Hel, and HellA4686 profiles to es- 
timate tblr, and then comparing the virial Mbh and the Mbh obtained assuming that the 
variable-fraction shifts observed in Mkn 1 10 are due to gravitational redshift, a measure- 
ment of the spin axis of the central black hole is possible 11121 : 6 ^ 21 it 5°. Unfortunately, 
there is as yet no convincing general evidence that the redward asymmetry seen in H/3bc 
profiles is entirely due to gravitational redshift ( ^4.4. II) . 

6 The Mass - Luminosity Diagram 

FigE] shows a plot of Mb vs. Mbh for a combined sample of 280 AGNs which cover 
an absolute B magnitude range -20 > Mb ^ - 27 and an Mbh interval 10^ < Mbh 
< 10^° (1 1431; cf. Refs. |281 41 1). The wide majority of sources are located between 
0.02 ^ L/LEdd ^ 1.00. RQ sources show evidence for significant Malmquist bias while 
RL sources show the opposite trend probably related to a bias towards selecting higher lu- 
minosity core-dominated sources likely to be beamed. We find that the Eddington limit 
defines an approximate upper boundary to the luminosity distribution |281 1, indicating that 
there are no \ow-z AGNs radiating significantly above the Eddington limit (this depends 
somewhat on the adopted Hq value). Figure also suggests that there may be fewer high 
(than low) mass sources with Lboi/MnH close to the Eddington limit. This might be a se- 
lection effect or an indication that galaxies with large Mbh may be unable to supply fuel 
at high Lijoi/Mbr i.e., a limit on M 11651 . If a sample comprising more high luminosity 
objects (L ^ 10^^ ergs s~^) is considered, 41% of the objects in that sample with L ^ 10^^ 
ergs s^^ have super-Eddington ratios |266|. We caution however that correction for ori- 
entation and beaming should decrease i^boi/^Edd- The lower envelope of the luminosity 
distribution of Fig. ^Jmay be due to a selection effect I281L or it may indicate that only 
sources radiating at 0.01 ^ (-Z^boi/^Edd) ~ 1-00 exhibit a stable BLR. The lower limit may 
be connected with the absence or small size of a standard accretion disk ( ^7.11 see also 
Refs. L113. 1651). Redshift evolution of the black hole mass-luminosity ratio is found to be 
much less than the (1 + z)^ evolution seen in QSO luminosity function evolution |48|. This 
means that quasars radiate in almost the same Lboi/MBH range also at high-z. 

7 The Eddington Ratio as the Driver of the El Correlations 

A PCA reveals a first principal component accounting for 48% of the spectrum-to-spectrum 
variance with just UV emission line strengths, ratios and widths as input. This linear combi- 
nation of the input variables depends on -Lboi/^Edd 121311279)1 . A reanalysis of PG quasars 
confirms that El is driven predominantly by LboiZ-^^Edd while a second principal component 
is driven by accretion rate 1211 . An UV based spectral PCA that uses the whole spectrum 
in the rest wavelength interval 1170-2100 A find consistent results although the relative 
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importance of luminosity and -Lboi/-^Edd is reversed i.e., iboi/-^Edd appears in a second 
eigenvector accounting for 20% of the spectrum-to-spectrum variance |285|. The El opti- 
cal sequence and the -^^boi/^Edd trends have been recently confirmed also for intermediate 
z quasars 11621 12331 128 71. Comparison between high and low-z QSOs confirms that the 
inverse Fellopt - [Olll]AA4959,5007 relationship (i.e., El) is indeed related to Lboi/^Edd 
rather than Mbh ll287l . 

Different Lboi/^BH values explain the change in Fell and W(H/3bc) prominence in 
the optical plane 11391 12781 11621 . The virial relationship can be rewritten in the form of 
FWHM(H/3bc) oc (Lboi/MBH)"^MBH"° '^ using the tblr- vLv correlation. The diag- 
nostic ratios Silll]A1892/ Clll]A1909 and AI111AI8O8/ Clll]A1909 provide an estimate of 
ne. Since correlates with FWHM(H/?bc) [278l ll39ll it is possible to write also the ion- 
ization parameter (and -RFeii) as a function of Lboi/MBH and Mbh 112891 . The El "elbow" 
sequence in the optical plane is reproduced by varying -Lboi/^Edd> with Lboi/-^^Edd — ^ 1 
toward spectral type A3. Orientation can be modelled as a third parameter in this scheme 
from the location of CD and LD RL QSOs I234II . Model computations suggest that sources 
of spectral type Al are actually almost face-on, low -Lboi/^Edd sources I139II . This result 
explains why some extremely low W(H/3bc) objects (the second extreme type mentioned 
in ^1.2t are observed in bin Al. Sources of bin Bl, A2, A3 are instead rather homogenous 
in terms of Lboi/^Edd with A/bh and Q acting as source od scatter. It has also been possi- 
ble to show that the H/3bc different profile shapes of Pop. A and B sources are linked to 
differences in Lboi/-Z^Edd II143I . 

The soft X-ray spectral index correlates extremely closely with Eddington ratio and 
FWHM(H/3bc) 1 88, , 229. .279 1. W(ClvA1549) decreases toward extreme population A (see 
a more detailed discussion in ^8.41 and Fig. where a dependence of W(ClvA1549) on 
spectral type is clearly shown). Also, the distributions of large CivA1549bc blueshifts and 
blue outliers seem to be governed by Lboi/^Edd ( lll43llT]1 '2881: Fig. flU i. Among type 1 
AGNs with large blueshifted [Olll]A5007, there is no correlation between the Eddington 
ratio and the amount of [Olll]A5007 blueshifts. However, the Eddington ratios of the blue 
outliers are the highest among AGNs with the same Mbh in a sample of 300 objects 1 143 1. 
These facts suggest that high Lboi/-^Edd is a necessary condition for [Olll]AA4959,5007 
large blueshifts. 

Both NLSyls galaxies and BAL QSOs of the PG sample he at the high Lboi/MeH 
extreme 1211 11391 1287L although these two object classes are well separated in a second 
eigenvector based on M I2TII . Other measurements confirm this scenario. BAL quasars 
show large ClvA1549 emission blueshifts 1 197','225V Using the luminosity and H/3bc to 
derive Mbh and m, both BAL and non-BAL QSOs at 2: ~ 2 tend to have even higher 
^boi/-^Edd than those at low z Extreme Pop. A sources may represent the most 

intense accretors with extreme BAL QSOs being observed at the largest inclination angles 
possible for a given wind geometry. A large mass flow may limit the viewing angle range 
of extreme Pop. A sources to be w 90° (S I3.1.1I & Fig. |9ll. It is important to stress a possible 
difference between the most extreme BALs (i.e., those with the highest terminal radial 
velocity, ~ 30000 km s~^) and the generality of BAL sources. While the most extreme 



162 



Paola Marziani, Deborah Dultzin-Hacyan and Jack W. Sulentic 



BAL QSOs seem to be associated to highest -Lboi/-^Edd> the presence of a shallow BAL 
of moderate terminal velocity is not necessarily a signature of high Lboi/-^^Edd since the 
terminal velocity correlates with Lboi/^Edd if the wind giving rise to the BALs is radiation 
driven l225lll96l fT97l. 

The separation between Pop. A and Pop. B sources, which may occur at Lboi/^Edd 
0.15 (Fig. fT3l could be in part explained by a change in accretion mode. We note 
that there is no correlation between radio loudness Rx and i^boi/-^Edd in a RL sample 
||92l although RL sources are differentiated in terms of Lboi/MBH from RQ sources 1 143 1. 
Most or all of the weakly AGNs in nearby galaxies are RL, highly sub-Eddington systems 
that are plausibly experiencing advection-dominated accretion [99 1. The vast majority of 
these objects appear to be radiating between 0.01 and 0.1 Eddington luminosity, and may 
therefore be akin to Pop. B powerful quasars. Radiatively inefficient accretion at less than 
a few per cent of the Eddington rate is unlikely to produce excess soft X ray emission 
11521 [159.1 . Lower mass loss rate, and the harder spectrum (which may decrease radiative 
acceleration) may both contribute to the absence of a strong wind in Pop. B sources. 

7.1 Are "Double Peakers" a Peculiar Population? 

Recent work seems to confirm the accretion disk origin of the LILs for double-peakers. 
Inter-profile comparison between Balmer lines and ClvA1549 (which, in a few cases, show 
a markedly different, single peaked Gaussian profile |68 1) can be interpreted in the frame- 
work of the disk and wind paradigm. More ClvA1549 - Balmer line profile intercompar- 
isons are needed to reach a firm conclusion also because the ClvA1549 profiles thus far 
observed cannot be immediately ascribed to a wind. The LIL emitting portion of the disk 
to produce double-peaks as separated as in Arp 102B is relatively modest, from ~ 100 i?g 
to ~ 500 Rg 1381 . At smaller radii, a hot solution may be appropriate 11591 . while at 
500 Rg self gravity may dominate. Using the customary empirical relation between the 
broad-line region size and optical continuum luminosity, Mbh and accretion rate have been 
computed for 135 AGNs with double-peaked broad emission lines from the SDSS from a 
survey of RL AGNs 11283 1 . The inferred range of Mbh of double-peakers cover almost 
the full range of Mbh of type-1 AGNs, 3-10^ Mq to 5-10^ Mgand Lboi up to 10^^ ergs 
s~^i283|. Typical Eddington ratios are 0.001, much lower than the typical type-1 AGN; 
about 90% of them show Lboi/^Edd ~ 0.02 (from the distribution of Lboi/^Edd shown by 
II283I '). Double-peaked sources lie close to the lower Ly^oi/MBu boundary or below in the 
mass-luminosity diagram. 

The evidence reported does not suggest that double peakers are truly peculiar sources, 
i.e., an independent Population C even if they lie out of the mass-luminosity sequence of 
Fig. A two-component model assuming that the line wings originate in an accretion 
disk with additional contribution from spherical emission region fits very well the observed 
H/3bc profiles of Pop. B sources fl87' '216 1. An evolutionary scheme from single peaked 
Pop. B and double peakers is at least conceivable: double peakers show, on average, the 
lowest i^boi/-^Edd (and probably also the lowest M). They may be AGNs exhausting their 
fuel supply, with little gas reservoir, in which only the innermost part of the accretion disk 
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Figure 14: Position of blue outliers (large blue spots) relative to other sources in our sample 
(filled circles). The largest circle indicates the position of PKS 0736+01, the only known RL 
blue outlier. The dashed line indicates the Eddington limit. Arrows indicate displacements 
of blue outliers if we apply an orientation correction 0.4 to the derived masses L140J . 

is left to produce broad emission lines. 

"Naked" AGNs Indirect arguments suggest that some "true" type 2 AGNs do exist 11271 
11651 . They are AGNs without a hidden BLR even in spectropolarimetric data |240|. The 
observed radius-luminosity relation for the BLR implies an increasing line width with de- 
creasing luminosity for a given black A/bh- However, the upper limit to the observed width 
of broad emission lines in AGNs at ~ 25000 km l229l may reflect a physical limit 
above which the BLR may not be able to survive 11271 . The disappearance of the BLR has 
been linked to the critical radius at which the disk changes from gas pressure dominated to 
radiation pressure dominated, or to disk evaporation 1 56. .164.1 . For low enough accretion 
rates, the critical radius becomes smaller than the innermost stable orbit of the accretion 
disk and the BLR may not form II165I . The intrinsic (i.e., unabsorbed) X-ray luminosity 
and an estimation of Mbh using the relationship between nuclear mass and bulge luminos- 
ity in galaxies indicate that all hidden BLR sources have accretion rates m >, 10^'^. This 
is also approximately the limit at which double-peakers are found L283J . Naked AGNs may 
this be observed if m<^ 10^^. 

7.2 On the Origin of Radio Loudness 

RL activity in low Mbh sources is observed, and there are RL NLSyls that radiate al 
logL/M « 4.5 II143I . There are no appreciable effects (within the limits set by S/N) on 
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the H/3bc profile attributable to radio loudness. Pop. B RL and RQ sources can have the 
same Mbh and Lboi/AfBH values. Therefore, it is not unreasonable to conclude that a 
similar range of Mbr and Lboi/AfeH i^ physically possible for both RQ and RL sources 
(2821. However, this does not mean that the mass function and the Lboi/AfBH distribution 
are necessarily the same for RL and RQ sources. RL and RQ sources are well separated in 
terms of Mbu in the El space 11211 at low z. A robust inference from a bootstrap analysis 
is that the mass function and the conditional probability of having certain Lboi/MeH val- 
ues at fixed Mbh are different for the two AGN classes fT43l. With the customary Mbh 
computation assumptions, Mbh in RL sources have been found to be Mbh ~ 10*^ M0, 
while only a few quasars may contain smaller black holes |92|. These conclusions are in 
agreement with studies based on the SDSS carried out using a sample of more than 6000 
AGNs [151]. RL sources were found to harbor systematically more massive black holes 
than are the RQ quasars with very high significance (see also Ref. |63 1). It is important that 
RQ and RL samples have indistinguishable distributions on the redshift-optical luminosity 
plane, excluding the possibility that either parameter is responsible for the observed Mbh 
difference 11431 11501 . Quasars from the FIRST Bright Quasar Survey 1 13] fill in the gap 
between the RL and RQ quasars in the radio vs. optical luminosity plane, show a contin- 
uous variation of radio luminosity with Mbh and no evidence for a discontinuity that may 
signal the turning on of powerful radio jets I124II . We suggest that several factors drive the 
RQ/RL separation in the El space: intrinsic mass function and Lboi/MBH distribution dif- 
ferences, jet-related effects yielding [Olll]AA4959,5007 enhancement in RL sources, and 
sample selection criteria 1 125, 21 143J. The range in radio luminosity at a given Mbh is 
several orders of magnitude. Additional parameters other than Mbh and Lboi/MBH like 
evolution should be invoked to explain the quasar radio-loudness dichotomy I15II . 

8 At the Origin and at the End 

8.1 Is the Emission Line Width Increasing with Luminosity? 

A trend between FWHM(H/?bc) and luminosity is expected if the usual virial and tblr- 
uL,y reasonings are applied. There is a well-defined lower boundary in the diagram Mbh 
vs. Mb (see Fig. |233|). Here we make the assumption that low-redshift NLSyl with 
the narrowest lines radiate very close to the Eddington ratio. If we assume logLboi/MBH 
PS 4.5, we obtain FWHM(H/3Bc)min ~ 600 km s^^ for logL = 11 in solar units and 
FWHMmin(H/3Bc)oc io(-oo8*^b) If vve consider the luminosity dependence of 
FWHM(H/3Bc)min> wc see that the expected trend for a = 0.6 reproduces fairly well 
the FWHM(H/3bc) lower boundary as a function of Mb- A less pronounced trend, es- 
pecially at high luminosity, is expected for a = 0.7. The presence of a correlation between 
FWHM(H/3bc) and luminosity in a survey may depend on sample and intrinsic dispersion 
of FWHM values in a narrow Mb range. One should consider that the expected luminosity 
dependence is very weak: AFWHM(H/3bc) increases by w 1000 km s^^ over an increase 
of AMb « 10, with FWHM(H^Bc)min changing from 1000 km s"^ to 2000 km s"^ In a 
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Figure 15: FWHM(H/3bc) in km s^ vs. absolute magnitude for two joined samples, 
one at low redshift (z 0.8, 11401 ). and the other for intermediate redshift quasars 12341 . 
Filled and open circles indicate RQ and RL sources respectively. The dependence on lu- 
minosity for the minimum FWHM(H/3bc) and for the boundary between Pop. A and B 
sources is shown. See text for details. 

narrow Mb range, the FWHM(H/3bc) intrinsic scatter spans the range 1000-10000 km s~^, 
therefore making any correlation intrinsically very weak even if statistically significant. A 
very large sample from the 2dF and 6dF redshift surveys shows a weak increase of H/3 line 
width with luminosity with a slope 0.15-0.2 (Av vs. log L), very close to the one derived 
from our calculation l48l . 

8.2 Lack of Evolution in the Fell Spectrum? 

An analysis of the spectra of 12 high-redshift quasars in the region of MgllX2800 shows 
remarkable similarity with the LBQS composite ^ 0.8; [237]). Near-infrared spectra of 
four QSOs located at z > 6 have been recently analyzed LI 03 J . Two out of four z > 6 QSOs 
have significant Felluv- Spectra of three additional sources indicate Felluv/MgllA2800 
similar to that observed in low-z quasars f78l l9^. The Felluv/MgllA2800 ratio at high 
redshift should be a factor of 3 lower than for low -redshift QSOs if the age of the universe 
at the earlier epoch is much less than 1 Gyr [,961. This prediction is due to the delayed 
contribution of type la supernovae to the iron abundance. Apart from the intriguing results 
for the highest z quasars, systematic examinations of the Felluv/MgllA2800 intensity ratio 
with redshift suggest that the ratio may indeed decline at z > 1.5 |104"1031. LBQS spectra 
show a significant Baldwin effect in Felluv emission. Further analysis reveals that the 
primary correlation of iron emission strength is probably with z, implying an evolutionary 
effect I^J. A very recent SDSS study based on ~ 16000 spectra in the reshift range from 
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0.08 to 5.41 shows that the dominant redshift effect is a result of the evolution of the blended 
Fell emission (optical) and the Balmer continuum 1285]. This redshift dependence can 
be explained by the evolution of chemical abundance in the quasar environment, or by 
an intrinsic change in the continuum itself |285|. In another study f25ri quasars were 
grouped on the basis of their Felluv/MgllA2800 intensity ratios. The fraction of quasars 
in each group with strong Felluv/MgllA2800 dropped significantly at z ^ 1.5. To check 
whether there is a change in iron abundance at high redshifts observations covering Fellopt 
are essential. 

8.3 A Scheme of Physical Evolution for the AGNs 

High-z quasars must have accreted at a relatively high rate in the early stages of their evo- 
lution, if black holes of Mbr ^ 10^ M© were already radiating at z ^ 6.4 12 57.26.6. 274J . 
From the low-z data, it is tempting to understand the Eddington ratio as the age of an AGN 
within its total lifetime (i.e., considering the totality of the time spent by a black hole as a 
"switched on" AGN): AGNs with steep X-ray spectra, strong Fell, and weak [Olll]A5007 
are AGNs in an early phase of their evolution, which proceeds roughly from the lower left 
corner (-Lboi/^Edd —>■ 1) of the optical plane of El toward the upper right corner (low 
Lboi/-A^BH> large Mbh)- In this hypothetical scenario local NLSyls are the "seedlings" 
of AGNs 12291 l88ll . Several results and correlations associated to the optical El can be 
explained by an evolutionary sequence. RL seem to be systematically more massive than 
RQ sources. They also radiate at lower Lboi/-^^Edd- They show extended NLRs with large 
W([Olll]AA4959,5007). At the other end we have "blue outhers" which are Pop. A sources 
with strong winds. Young AGNs still possess a compact NLR; they may have plenty of fuel 
and may show spectral evidence of circum-nuclear Starbursts. More massive, evolved sys- 
tems may have less fuel and have had all the time to ionize interstellar clouds in the bulge 
of their host galaxies as well as to produce radio lobes exceeding the host galaxy size. Dy- 
ing quasars have a limited gas supply surrounding their continuum-emitting regions, and 
may at best produce emission lines through the illumination of a last "strip" or ring of op- 
tically thick material as it may be the case of double-peakers. In this sense, the iboi/-^Edd 
sequence of El becomes literally an age sequence. 

8.4 The Baldwin Effect 

It is important to stress that the Baldwin effect is a very loose correlation between specific 
luminosity and ClvA1549 equivalent width. Claims and counter-claims of a Baldwin ef- 
fect on the basis of small samples (few tens of objects) are unreliable 1229 1; the statistical 
weakness of the Baldwin correlation implies that the effect becomes appreciable only if a 
very large range in luminosity is considered, 4-6 decades. As pointed out by Sulentic et 
al. f2291, results until mid- 1999 have led to a standard scenario in which the Baldwin effect 
occurs in all measurable HILs except NvA1240, and the slope of the Baldwin relationship 
(i.e., W(ClvA1549)oc L~'^, |5|) increases with ionization potential (e.g. |168|; see Ref. 
l96t for the implications of the apparent lack of any Baldwin effect for NvA1240). These 
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Figure 16: The Baldwin effects in the HST based sample of Bachev et al. EJ. Left panel: 
dependence of rest-frame W(CivA1549bc) on bolometric luminosity. Middle panel: de- 
pendence of W(civbc) on the L/M ratio (in solar units). Filled circles indicate RQ, open 
circles RL sources. Right Panel, same as middle panel but with median averages with the 
spectral types defined according to Ref. |230|. 

results have been basically confirmed by more recent studies based on large quasar samples 
||54l|59l[l98|, LBQS |87| and SDSS |285| included. A source of dispersion involves low- 
luminosity objects that have low W(ClvA1549) ~ 10 30 A since they tend to blur even 
more the Baldwin relationship. These sources are NLSyls that preferentially occupy the 
lower left part of a log W vs. log L diagram. 

Claims of a Baldwin effect on MgllA2800 should be taken somewhat with care f2851, 
since the underlying Felluv emission can be very strong as well as variable from source to 
source depending on the BLR physical conditions, and the Felluv behavior has a function 
of luminosity and redshift is still poorly understood. Similar claims for W(H/3bc) have 
been rather erratic and unconvincing due to small sample sizes and large intrinsic scatter 
(Refs. 1 146. 229] found no trend; Ref. 1 162 273 1 did). 

The [Olll]AA4959,5007 lines seem to decrease in equivalent width with increasing lu- 
minosity I59lll62ll233ll285]l . In the joint data of Refs. l233ll and fl40ll there is definitely a 
significant correlation. Even if a luminosity correlation may not be always detectable (we 
stress again that the luminosity range must be very large) weak or no [Olll]AA4959,5007 
seems to be a common property of many luminous AGNs 11621 12871 . The reason why 
[Olll]AA4959,5007 seem to experience a Baldwin effect is presently not clear. It is possi- 
ble that the [Olll] AA4959,5007 luminosity is limited by the physical size of the NLR |1621, 
or that we are seeing an evolutionary effect at high-z, making high-z sources more similar 
to low-z Pop. A AGNs. 

8.4.1 The Baldwin Effect as an Evolutionary Effect 

A recent analysis of 80 PG quasars and one of w 120 HST archive spectra indicates a strong 
correlation of W(ClvA1549) with some of the emission parameters which define El. Since 
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^boi/^Edd drives the El correlations, Lboi/^Edd may be the primary physical parameter 
behind the Baldwin effect for ClvA1549 (|4 flTlEll; Fig. [T6ll. Also, high-z quasars show 
large ClvA1549 blueshifts like Pop. A sources at low-^ lT98ll233ll225l . Near-IR, low res- 
olution spectra of eight of the most distant quasars (in the range 4.9 < z < 6.4) show that 
half of these quasars are characterized by deep, broad and blueshifted absorption features 
typical of BAL quasars II135I . Similarities and close association between NLSyl, BAL 
QSOs, strong Fell emission, and low W(ClvA1549) led to the suggestions that NLSyls 
might be Seyfert galaxies in their early stage of evolution and as such may be low-redshift, 
low-luminosity analogues of high-redshift quasars 1 144, 229|. The ClvA1549 profiles show 
decreasing equivalent widths with increasing luminosity and the line profiles become sim- 
ilar to the "trapezoidal" shapes often observed in low-z NLSyls I266II . However, a direct 
comparison of composite spectra of NLSyl and z ~ 4 quasars through a PC A suggests 
that high-z quasi-stellar objects do not show a strong preference toward NLSyl behavior 
lEr6 |. Since this result is sample dependent and the mass-luminosity diagram of quasars can 
be similar up to z ~ 2.5 ll48l . it is still debatable whether evolution of Lboi/-^^Edd with z, 
selection effects or a combination of both may lead to the Baldwin effect [2291. After all, 
we are not observing other evolutionary effects even at z w 6. What seems unlikely is that 
the driving parameter of the Baldwin effect is Mbh 12651 . Computing A/bh employing 
FWHM(ClvA1549) and assuming virial motions does not demonstrate that Mbh is driving 
the Baldwin effect because it is assumed to be Mbhoc L^'^ , and this assumption leads, in 
part, to circular reasoning. 

8.5 Quasars as Standard Candles 

The Supernova Cosmology Project results supporting a cosmological constant A 7^ are 
still highly uncertain. The expected effect is just a few tens of magnitude at z w 0.6, and 
goes down to at z 1.5. A robust verification of a cosmological model with relative en- 
ergy density JIm ~ 0.28 and I^a ~ 0.72 would require accurate measurements of standard 
candles in the range 0.6 ^ z 3, as it can be seen from the Hubble diagram with type la 
supernovae 1200., 201 1. 

It is clear from the previous discussion that, at high z, we are observing quasars that 
can be very similar to the AGNs we are observing at low z, in terms of line width, Fell 
prominence, HIL equivalent widths. Trends with Lboi/MBH are much better defined than 
trends with luminosity. Luminosity effects remain weak and prone to sample bias (see, for 
example, Fig. where the Baldwin effect is not significant while W(ClvA1549) shows 
a significant correlation with -f/boi/^Edd)- Furthermore, Lboi/-^^Edd seems to affect many 
parameters related to physics and structure of the BLR, from LIL widths to HIL shifts and 
ionization degree. It makes therefore sense to use El correlations to estimate Lboi/A/BH- If 
Mbh is known with reasonable accuracy, it may become possible to retrieve z-independent 
information on L^oi- In principle, this can be achieved through two methods: 

• a multivariate analysis which isolates the best linear combination of variables which 
are dependent only on Lboi/MBH and Mbh I279II : 
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• an extension of the El diagram to high luminosity and high z I141II . 

Instruments like the infrared spectrometers mounted on 10-m sized telescopes make 
possible observations of intermediate redshift quasars (1 ^ z ^ 2.2) with SjN and reso- 
lution sufficient to apply the same data analysis procedure used for optical spectra of sources 
with z If these data are available, it is possible to measure the El parameters and to 
build an extension of "optical" El diagram. To fully exploit the El diagram, orientation 
needs to be estimated on an object-by-object basis. In the case of a radiation pressure 
driven wind, the ClvA1549 shift and profile shape should be sensitive to both Lboi/MeH 
and Q. Therefore, ClvA1549 measurements may provide an independent estimate of ori- 
entation, and hence a 3D observational parameter space to map into a 3D physical space 
{Q, Lboi/-^BH> -^bh) which will ultimately yield a redshift-independent estimate of the 
luminosity (if an accurate low-z calibration can be defined). Merging results on Mbh and 
on Lboi from the previous analysis, it is reasonable to infer that, if Q is known within ±5°, 
errors could be ~ 60% on Mbh> 80% on Lboi/AfeH- If errors are distributed randomly, 
an uncertainty as high as Am 1 could be farther reduced if samples of sources within a 
spectral type are averaged, yielding valuable data points for the Hubble diagram. 



9 Conclusion 

The main purpose of this paper was to point out how it is becoming possible to measure 
accretion parameters (especially Mbh and Lboi/-^^BH) of AGNs. We highlighted major 
results and sources of uncertainty, as well as the connections between Mbh and -Lboi/^Edd, 
the structure of the broad line emitting region, and physical evolution. If we compare our 
present knowledge to the one of about 10 years ago, we can see an impressive progress. 
Mbh determinations through the virial assumption have become widely popular, and have 
been applied to quasars at all known z, up to 6.4, even if several sources of uncertainty limit 
strongly the accuracy by which Mbh is known for individual objects. It has been possible 
to show that many line properties and trends (including the Baldwin effect) correlate with 
Lboi/-^Edd> and the the BLR structure is strongly affected by LboiZ-i^Edd through El. There 
is no doubt that we have a much better framework with a strong observational support to 
understand the BLR structure, minority objects like "double-peakers" and BAL QSOs, and 
source evolution. Improvements on Mbh and Lboi/MnH estimates are possible through 
new multivariate analysis including orientation sensitive parameters, or through disk and 
wind BLR models accounting for the combined effects on line parameters of orientation and 
physics. We may be able to achieve a satisfactory vision of quasar structural evolution (not 
to mention boundaries on cosmography) well before the emitting regions may be resolved 
through space interferometry missions (e.g., Ref. [245 1). 

We thank Dr. Massimo Calvani for a careful reading of the manuscript. 
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